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^ ' We review a variety of theoretical and experimental investigations aimed at 

■ proving our knowledge of the nuclear matter equation of state. Of particular interest 
are nuclear matter extreme states in terms of density and/or isospin asymmetry. The 
equation of state of matter with unequal concentrations of protons and neutrons has nu- 
merous applications. These include heavy-ion collisions, the physics of rare, short-lived 
nuclei and, on a dramatically different scale, the physics of neutron stars. The "common 
denominator" among these (seemingly) very different systems is the symmetry energy, 
which plays a crucial role in both the formation of the neutron skin in neutron-rich 
nuclei and the radius of a neutron star (a system 18 orders of magnitude larger and 55 
orders of magnitudes heavier). The details of the density dependence of the symmetry 

■ energy are not yet sufficiently constrained. Throughout this article, our emphasis will 
f***^ ' be on the importance of adopting a microscopic approach to the many-body problem, 

' which we believe to be the one with true predictive power. 
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1. Introduction 

The interaction of hadrons in nuclei is a problem that goes to the very core of 
nuclear physics. In fact, our present knowledge of the nuclear force in free space is, 
in itself, the result of decades of struggle,^ which will not be reviewed in this article. 
The nature of the nuclear force in the medium is of course a much more complex 
problem, as it involves aspects of the force that cannot be constrained through 
free-space nucleon-nucleon (NN) scattering. Predictions of properties of nuclei are 
the ultimate test for many-body theories. 

Nuclear matter is an alternative and convenient theoretical laboratory for many- 
body theories. By "nuclear matter" we mean an infinite system of nucleons acted on 
by their mutual strong forces and no electromagnetic interactions. Nuclear matter 
is characterized by its energy /particle as a function of density and other thermody- 
namic quantities, if appropriate (e.g. temperature). Such relation is known as the 
nuclear matter equation of state (EoS). The translational invariance of the system 
facilitates theoretical calculations. At the same time, adopting what is known as 
"local density approximation" , one can use the EoS directly in calculations of fi- 
nite systems. This procedure is applied, for instance, in Thomas-Fermi calculations 
within the liquid drop model, where an appropriate energy functional is written in 
terms of the EoS.^EEl 

Symmetric nuclear matter (that is, matter with equal densities of protons and 
neutrons) has been studied extensively. The so-call ed conventional approach goes 
back to earlier works by Brueckner and 

others^™ 

and is known as the Brueckner- 
Hartree-Fock (BHF) theory. The Brueckner theory is based on a li nked-cluster per- 
turbation series of the ground state energy of a many-body system!S) 6 | 7 | 8|9|10 ] Such 
series was shown to converge when the cluster diagrams are regrouped according 
to the number o f hole lines. The variational approach was also pursued as an al- 
ternative 

methodULL^ 

and yielded predictions in close agree ment with those from 
Brueckner theory when realistic NN potentials were employed 

Durin g the 198 0's, the Dirac-Brueckner-Hartree-Fock (DBHF) approach was 
developedPHSllfil The main break-through came with the observation that the 
DBHF theory, unlike conventional Brueckner theory, could describe successfully 
the saturation properties of nuclear matter, that is, saturation energy and density 
of the equation of state. The DBHF method adopts realistic NN interactions and 
contains important relativistic features. It describes the nuclear mean field in terms 
of strong, competing scalar and vector fields that, together, account for the binding 
of nucleons as well as the large spin-orbit splitting seen in nuclear states. 

Regardless of the chosen many-body theory, a quantitative NN potential should 
be part of its input. Recently, chiral effective theories of the nuclear fore have be- 
come popular as a mean to respect the symmetries of Quantum Chromo Dynamics 
(QCD) while retaining the basic degrees of freedom typical of low-energy nuclear 
physics. Chiral effective theories provide a well-defined scheme to determine the 
appropriate many-body diagrams to be included at each order of the perturbation. 
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However, being based on a low-momentum expansion, interactions derived from chi- 
ral perturbation theory are not suitable for applications in dense nuclear/neutron 
matter, where high Fermi momenta are involved. Instead, meson-theoretic or phe- 
nomenological NN potentials are typically employed as input to the many-body 
theory. Mean-field models (relativistic and non-relativistic) are a popular, although 
non-microscopic, alternative to methods based on the in-medium reaction matrix, 
such as BHF and DBHF. 

It is the purpose of this article to review the status of microscopic studies of 
nuclear and neutron-rich matter, with particular emphasis on the latter, as it re- 
lates to present empirical investigations within terrestrial laboratories (heavy-ion 
collisions) or astrophysical observations (neutron stars). At the same time, we will 
review the present status concerning available empirical information which can be 
utilized to guide and constrain theories. We also wish to provide a self-contained 
account of the recent work with asymmetric matter done by the Idaho group, in- 
cluding extensive numerical tables for the interested user. 

Clearly, the goal to describe the properties of (dense) many-body systems con- 
sistently from the underlying forces and including all potentially important mech- 
anisms, is a most ambitious program and far from having been completed. The 
importance of pursuing a microscopic approach towards the accomplishment of 
this goal is a theme that will surface repeatedly throughout this article. 

2. Asymmetric Nuclear Matter: Some Facts and Phenomenology 

Isospin-asymmetric nuclear matter (IANM) simulates the interior of a "nucleus" 
with unequal densities of protons and neutrons. The equation of state of (cold) 
IANM is then a function of density as well as the relative concentrations of protons 
and neutrons. 

The recent and fast-growing interest in IANM stems from its close connection 
to neutron-rich nuclei, or, more generally, asymmetric nuclei, including the very 
"exotic" ones known as "halo" nuclei. At this time, the boundaries of the nuclear 
chart are uncertain, with a few hundreds stable nuclides known to exist and per- 
haps a few thousands believed to exist. The Facility for Rare Isotope Beams (FRIB) 
has recently been approved for design and construction at Michigan State Univer- 
sity (MSU). The facility will deliver intense beams of rare isotopes, the study of 
which can provide crucial information on short-lived elements normally not found 
on earth. 10 Thus, this new experimental program will have widespread impact, 
ranging from the origin of elements to the evolution of the cosmos. 

It is estimated that the design and construction of FRIB will take ten yearsP^In 
the meantime, systematic investigations to determine the properties of asymmetric 
nuclear matter are proliferating at existing facilities. ^^^^ 

From the theoretical side, some older studies of IANM can be found in RefsP222Q] 
Interactions adjusted to fit properties of finite nuclei, such as those based on the non- 
relativistic Skyrme Hartree-Fock theorjEH or the relativistic mean field theory)^ 
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have been used to extract phenomenologicalEoS. Variational calculations of asym- 
metric matter have also been reportedP^ES] Fuchs et a/P31 defined a Lorentz in- 
variant functional of the baryon field operators to project Dirac-Brueckner nuclear 
matter results onto the meson-nucleon vertices of an effective density-dependent 
field theory. This was then applied to asymmetric matter and finite nuclei in Hartree 
calculations E51 E xtensiv e work with IANM has also been reported by Lombardo 
and collaborators!2SI23 Dirac-Brueckner- Hartree- Fock calculations of IANM prop- 

IOQI I9QI 

erties were performed by the Oslo group, the Idaho group, 1 ^ 1 and by Fuchs and 
collaborators Typically, considerable model dependence is observed among the 
different EoS of IANM, especially in the high-density region. 

Asymmetric nuclear matter can be characterized by the neutron density, p n , 
and the proton density, p p , defined as the number of neutrons or protons per unit 
of volume. In infinite matter, they are obtained by summing the neutron or proton 
states per volume (up to their respective Fermi momenta, kp or k p F ) and applying 
the appropriate degeneracy factor. The result is 

Of) 3 (1) 
Pl ~ 3tt 2 ' 1 ' 

with i — n or p. 

It is more convenient to refer to the total density p = p n + p P and the asymmetry 
(or neutron excess) parameter a = p " - Pp . Clearly, a—0 corresponds to symmetric 
matter and a=l to neutron matter. In terms of a and the average Fermi momentum, 
&F, related to the total density in the usual way, 

the neutron and proton Fermi momenta can be expressed as 

k F = k F (l + a) 1/3 (3) 

and 

k p F = Mi - «) 1/3 , (4) 

respectively. 

The energy /particle in IANM can, to a very good degree of approximation, be 
written as 

e(p,a) w e (p) + e sym (p)a 2 , (5) 

where the first term is the energy /particle in symmetric matter and e sym is known 
as the symmetry energy. In the Bethe-Weizsacker formula for the nuclear binding 
energy, it represents the amount of binding a nucleus has to lose when the numbers 
of protons and neutrons are unequal. The symmetry energy is also closely related 
to the neutron /3-decay in dense matter, whose threshold depends on the proton 
fraction. A typical value for e sym at nuclear matter density (po) is 30 MeV, with 
theoretical predictions spreading approximately between 26 and 35 MeV. The effect 
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of a term of fourth order in the asymmetry parameter (C(a 4 )) on the bulk properties 
of neutron stars is very small, although it may impact the proton fraction at high 
density. More generally, non-quadratic terms are usually associated with isovector 
pairing, which is a surface effect and thus vanishes in infinite matterPH 

Equation (|5|) displays a convenient separation between the symmetric and aym- 
metric parts of the EoS, which facilitates the identification of observables that may 
be sensiti ve, for i nsta nce, mai nly to the symmetry energy. At this time, groups 



be sensiti ve, tor instance, mainly to the symmetry energy. At this time, groups 
from GSIpS MSUpS ItalyJSl FranceJSS and ChinEp^SZEl are investigating the 
density dependence of the symmetry energy through heavy-ion collisions. From re- 
cent results, such as those reported at the 2009 "International Workshop on Nuclear 
Dynamics in Heavy-ion Reactions and the Symmetry Energy" (Shanghai, China, 
August 22-25, 2009), these investigations appear to agree reasonably well on the 
following paramctrization of the symmetry energy: 

e sym {p) = 12.5Mev(^) 2/3 + l7.5Mev(^Y\ (6) 

v PO ' X PQ J 

where the first term is the kinetic contribution and ji (the exponent appearing in 
the potential energy part) is found to be between 0.4 and 1.0. Naturally, there are 
uncertainties associated with all transport models. Recent constraints from MSlJ^ 
were extracted from simulations of 112 Sn and 124 Sn collisions with an Improved 
Quantum Molecular Dynamics transport model and are consistent with isospin 
diffusion data and the ratio of neutron and proton spectra. 

Typically, paramctrizations like the one given in Eq. ([6]) are valid at or below 
the saturation density, po- Efforts to constrain the behavior of the symmetry energy 
at higher densities are presently being pursued through observables such as 7r~ /tt + 
ratio, K + /K° ratio, neutron/proton differential transverse flow, or nucleon elliptic 
flowESI 

A more detailed discussion of the symmetry energy and the status of its theo- 
retical predictions will be presented in Sec. 3.2.3. 



3. The Microscopic Approach 

In this section, we present a discussion of the microscopic approach to nuclear 
matter, in general, and the DBHF method, in particular. 



3.1. The two-body potential 

By ab initio we mean that the starting point of the many-body calculation is a 
realistic NN interaction which is then applied in the nuclear medium without any 
additional free parameters. Thus the first question to be confronted concerns the 
choice of the "best" NN interaction. As already mentioned in the Introduction, after 
the development of QCD and the understanding of its symmetries, chiral effective 
theories^ were developed as a way to respect the symmetries of QCD while keeping 
the degrees of freedom (nucleons and pions) typical of low-energy nuclear physics. 
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However, chiral perturbation theory (ChPT) has definite limitations as far as the 
range of allowed momenta is concerned. For the purpose of applications in dense 
matter, where higher and higher momenta become involved with increasing Fermi 
momentum, NN potentials based on ChPT are unsuitable. 

Relativistic meson theory is an appropriate framework to deal with the high mo- 
menta encountered in dense matter. In particular, the one-boson-exchange (OBE) 
model has proven very successful in describing NN data in free space and has a 
good theoretical foundation. Among the ma ny availab le OBE potentials, some be- 
ing part of the "high-precision generation" | 40 | 41 | 42 [ we see k a mom entum-space 
potential developed within a relativistic scattering equation, such as the one ob- 
tained through the ThompsorPSl three-dimensional reduction of the Bethe-Salpeter 
equation!^ Furthermore, we require a potential that uses the pseudovector coupling 
for the interaction of nucleons with pseudoscalar mesons. With these constraints in 
mind, as well as the requirement of a good description of the NN data, Bonn is a 
reasonable choice. As is well known, the NN potential model dependence of nuclear 
matter predictions is not negligible. The saturation points obtained with different 
NN potentials move along the famous "Coester band" depending on the strength 
of the tensor force, with the weakest tensor force yielding the largest attraction. 
This can be understood in terms of medium effects (particularly Pauli blocking) 
reducing the (attractive) second-order term in the expansion of the reaction matrix. 
A large second-order term will undergo a large reduction in the medium. Therefore, 
noticing that the second-order term is dominated by the tensor component of the 
force, nuclear potentials with a strong tensor component will yield less attraction 
in the medium. For the same reason (that is, the role of the tensor force in nuclear 
matter), the potential model dependence is strongly reduced in pure (or nearly 
pure) neutron matter, due to the absence of isospin-zero partial waves. 

Already when QCD (and its symmetries) were unknown, it was observed that 
the contribution from the nucleon-antinucleon pair diagram, Fig. [TJ becomes un- 
reasonably large if the pseudoscalar (ps) coupling is used, leading to very large 
pion-nucleon scattering lengths!^ We recall that the Lagrangian density for pseu- 
doscalar coupling of the nucleon field (ip) with the pseudoscalar meson field (<f>) 
is 

C ps = -ig pa ipj 5 ip(f). (7) 

On the other hand, the same contribution (Fig. [IJ is heavily reduced by the pseu- 
dovector (pv) coupling (a mechanism which became known as "pair suppression" ) . 
The reason for the suppression is the presence of the covariant derivative at the 
pseudovector vertex, 

C pv = ^^757 M W, (8) 

which reduces the contribution of the vertex for low momenta and , thus, explains 
the small value of the pion-nucleon scattering length at threshold!^ Considerations 
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Fig. 1. Contribution to the NN interaction from virtual pair excitation. Upward- and downward- 
pointing arrows represent nucleons and antinucleons, respectively. Dashed lines denote mesons. 



based on chiral symmetry^ can further motivate the choice of the pseudovector 
coupling. 

In closing this section, we wish to highlight the most important aspect of the 
"ab initio" approach: namely, the only free parameters of the model (the param- 
eters of the NN potential) are determined by fitting the free-space NN data and 
never readjusted in the medium. In other words, the model parameters are tightly 
constrained and the calculation in the medium is parameter free. The presence of 
free parameters in the medium would generate effects and sensitivities which are 
hard to control and interfere with the predictive power of the theory. 



3.2. The Dirac-Brueckner-Hartree-Fock approach to symmetric 
and asymmetric nuclear matter 

3.2.1. Formalism 

The main strength of the DBHF approach is in its inherent ability to account for 
important three-body forces through its density dependence. In Fig. [2] we show a 
three-body force (TBF) originating from virtual excitation of a nucleon-antinucleon 
pair, known as "Z-diagram" . Notice that the observations from the previous section 
ensure that the corresponding diagram at the two-body level, Fig. [TJ is moderate 
in size when the pv coupling is used. The main feature of the DBHF method turns 
out to be closely related to the TBF depicted in Fig. [21 as we will argue next. In the 
DBHF approach, one describes the positive energy solutions of the Dirac equation 
in the medium as 



,'E* +m*\ 1/2 / 1 \ 



where the effective mass, m* , is defined as m* = m+Us, with Us an attractive scalar 
potential. (This will be derived below.) It can be shown that both the description of 
a single-nucleon via Eq. ([9]) and the evaluation of the Z-diagram, Fig. [21 generate a 
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Fig. 2. Three-body force due to virtual pair excitation. Conventions as in the previous figure. 

repulsive effect on the energy/particle in symmetric nuclear matter which depends 
on the density approximately as 

, (10) 

and provides the saturating mechanism missing from conventional Brucckner cal- 
culations. (Alternatively, explicit TBF are used along with the BHF method in 
order to achieve a similar result.) Brown showed that the bulk of the desired effect 
can be obtained as a lowest order (in p 2 /m) relativistic correction to the single- 
particle propagation.SSl With the in-medium spinor as in Eq. ([9]) , the correction to 
the free-space spinor can be written approximately as 

«*(p,A)-ti(p,A)« ( a % ) xx, (11) 

where for simplicity the spinor normalization factor has been set equal to 1, in 
which case it is clearly seen that the entire effect originates from the modification 
of the spinor's lower component. By expanding the single-particle energy to or- 
der {/§, Brown showed that the correction to the energy consistent with Eq. 
can be written as ^(^f) 2 - He then proceeded to estimate the correction to the 
energy /particle and found it to be approximately as given in Eq. (fTDj) . 

The approximate equivalence of the effective-mass description of Dirac states 
and the contribution from the Z-diagram has a simple intuitive explanation in the 
observation that Eq. like any other solution of the Dirac equation, can be writ- 
ten as a superposition of positive and negative energy solutions. On the other hand, 
the "nucleon" in the middle of the Z-diagram, Fig. [3J is precisely a superposition 
of positive and negative energy states. In summary, the DBHF method effectively 
takes into account a particular class of TBF, which are crucial for nuclear matter 
saturation. 

Of course, other, more popular, three-body forces (not included in DBHF) need 
to be addressed as well. Figure [3] shows the TBF that is included in essentially all 
TBF models, regardless other components; it is the Fujita-Miyazawa TBF.S3 With 
the addition of contributions from ttN S-waves, one ends up with the well-known 
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Fig. 3. Left: three-body force arising from A-isobar excitation (thick line). Right: two-meson 
exchange contribution to the NN interaction involving A-isobar excitation. 



Tucson-Melbourne TBF.SS1 The microscopic TBF of RefPEl 

includes contributions 

from excitation of the Roper resonance (Pn isobar) as well. 

Now, if diagrams such as the one shown on the left-hand side of Fig. |3] are 
taken into account, consistency requires that medium modifications at the corre- 
sponding two-body level are also included, that is, the diagram on the right-hand 
side of Fig. [3] should be present and properly medium modified. L arg e cancella- 
tions then take place, a fact that was brought up a long time ago!^ When the 
two-body sector is handled via OBE diagrams, the two-pion exchange is effectively 
incorporated through the a "meson" , which cannot generate the (large) medium 
effects (dispersion and Pauli blocking on A intermediate states) required by the 
consistency arguments presented above. Thus, caution needs to be exercised when 
applying TBF in a particular space (of nucleons only, or nucleons and A). 

Having first summarized the mai n DBHF philosophy, we now proceed to describe 
the DBHF calculation of IANm!^ In the end, this will take us back to the crucial 
point of the DBHF approximation, Eq. (J9)). 

We start from the Thompson ^ relativistic three-dimensional reduction of the 
Bethe-Salpeter equation!^ The Thompson equation is applied to nuclear matter in 
strict analogy to free-space scattering and reads, in the nuclear matter rest frame, 

9ij((f, q, P, (e*j)o) = v* 3 {q, q) 

,i2) 

where gij is the in-medium reaction matrix (ij=nn, pp, or np), and the asterix 
signifies that medium effects are applied to those quantities. Thus the NN potential, 
v*j , is constructed in terms of effective Dirac states (in- medium spinors) as explained 
above. In Eq. (1121) . q, <f , and K are the initial, final, and intermediate relative 
momenta, and E* = J (m*) 2 + K 2 . The momenta of the two interacting particles 
in the nuclear matter rest frame have been expressed in terms of their relative 
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momentum and the center-of-mass momentum, P, through 

P = k l + k 2 , (13) 

K = ^. (14) 

The energy of the two-particle system is 

e* j (P,K) = e*(P,K) + e*(P,K) (15) 

and (e^o is the starting energy. The single-particle energy e* includes kinetic en- 
ergy and potential energy contributions (see Eq. (|2"91 below). The Pauli operator, 
Qij, prevents scattering to occupied nn, pp, or np states. To eliminate the angular 
dependence from the kernel of Eq. (TT21 , it is customary to replace the exact Pauli 
operator with its angle- aver age. Detailed expressions for the Pauli operator and the 
average cent er-o f-mass momentum in the case of two different Fermi seas can be 
found in RefP. 



With the definitions 



g« = t3^tSW (is) 



and 



m 3 



V* = —v* J — , (17) 



one can rewrite Eq. (|12l) as 

G^,q,P, {e* j ) Q ) = V*0,q) 



W W> g) (^-V.O Ctf(Klg (eli,0) • (18) 

which is formally identical to its non-relativistic counterpart. 

The goal is to determine self-consistently the nuclear matter single-particle po- 
tential which, for IANM, will be different for neutrons and protons. To facilitate 
the description of the procedure, we will use a schematic notation for the neu- 
tron/proton potential. We write, for neutrons, 

and for protons 

U p = Up„ + U pp , (20) 

where each of the four pieces on the right-hand-side of Eqs. (|19H20p signifies an 
integral of the appropriate G- matrix (nn, pp, or np) obtained from Eq. (fT8|) . Clearly, 
the two equations above are coupled through the np component and so they must 
be solved simultaneously. Furthermore, the G-matrix equation and Eqs. (|19H20[) 
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are coupled through the single-particle energy (which includes the single-particle 
potential, itself defined in terms of the G-matrix) . So we have a coupled system to 
be solved self-consistently. 

Before proceeding with the self-consistency, one needs an ansatz for the single- 
particle potential. The latter is suggested by the most general structure of the 
nucleon self-energy operator consistent with all symmetry requirements. That is: 

Ui(f) = UsAp) + 7oE#» - 7 • pUvAp) > (21) 

where Us,i and Uy,i are an attractive scalar held and a repulsive vector field, re- 
spectively, with U Vi the timelike component of the vector field. These fields are in 
general density and momentum dependent. We take 

U t (p)~UsAp)+loU v ^(p), (22) 

which amounts to assuming that the spacelike component of the vector field is 
much smaller than both Us,i and E/y-j. Furthermore, neglecting the momentum 
dependence of the scalar and vector fields and inserting Eq. (|22|) in the Dirac 
equation for neutrons/protons propagating in nuclear matter, 

(TmP* 1 -TTH- Ui(p))ui{fr A) = , (23) 

naturally leads to rewriting the Dirac equation in the form 

(7^*-m?H(p,A)=0, (24) 

with positive energy solutions as in Eq. m* = m + Us,%, and 

( P °y=p -u°Ap). (25) 

The subscript "i" signifies that these parameters are different for protons and neu- 
trons. 

As in the symmetric matter case evaluating the expectation value of Eq. (|22jl 
leads to a parametrization of the single particle potential for protons and neutrons 
(Eqs. (fT9ll20|) ) in terms of the constants Us,i and Uy i which is given by 

TYi * 777 * 

u^) = -=£< mm* >= -j£u s ,i + u Vti . (26) 

i i 

Also, 

Ui(p)= G iM,f?j)i (27) 

which, along with Eq. (f2l)|). allows the self-consistent determination of the single- 
particle potential as explained below. 

The kinetic contribution to the single-particle energy is 

Ti(p) = <p\"/-p + m\p>= ^r-^, (28) 
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and the single-particle energy is 

e* (p) = Ti (p) + Ui (p) = El + U$ ti . (29) 

The constants m* and 

U ,i = U s ,i + U$ ti (30) 

are convenient to work with as they facilitate the connection with the usual non- 
relativistic frameworkP 

Starting from some initial values of m* and J7o,i, the G-matrix equation is solved 
and a first approximation for Ui (p) is obtained by integrating the G-matrix over the 
appropriate Fermi sea, see Eq. (|27p. This solution is again parametrized in terms of 
a new set of constants, determined by fitting the parametrized Ui, Eq. (|26j) . to its 
values calculated at two momenta, a procedure known as the "reference spectrum 
approximation" . The iterative procedure is repeated until convergence is reached. 

Finally, the energy per neutron or proton in nuclear matter is calculated from 

Si = j < Ti > +^ <U l >-m. (31) 
The EoS, or energy per nucleon as a function of density, is then written as 

e{Pn,Pp) = — , (32) 

or 

, (1 + a)e n + (1 - a)e p 
e(k F ,a) = . (33) 

Clearly, symmetric nuclear matter is obtained as a by-product of the calculation 
described above by setting a=0. 

In the DBHF calculation of RefPQl, a similar scheme is applied to obtain the 
self-consistent G-matrix and spinor basis. At each step of the iterative procedure, 
the nucleon self-energy is first calculated by integrating the G-matrix elements over 
the Fermi sea. The self-energy components are then obtained from the appropriate 
traces^! i. e. 



E s = itr[£], E = ~tr[7b£], = -ttUM7 • ■ (34) 
4 4 i\p] z 

A covariant representation of the G-matrix is used to facilitate the transitions 
between the two-nucleon center-of-mass frame and the nuclear matter rest frame. 
Again, the pre sence of the medium naturally leads to define an effective mass, which 
is given by^ 

m*(p,k F ) = (m + ReX s (p,k F ))/(l + X v (p,k F )). (35) 

Although the effective mass is density and momentum dependent, ultimately it 
is defined as the value of the expression given in Eq. ([33]) at p = k F . Thus, the 
"reference spectrum approximation" is also employed in Ref|3ill, although differently 
in some of the technical aspects as compared to the Idaho method. These technical 
aspects appear to impact the predictions of some isovector quantities, as it will be 
discussed in Sec. 3.2.4. 
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0.5 r 1.5 

P (fm- 3 ) 

Fig. 4. DBHF predictions for the EoS of symmetric matter (solid red) and neutron matter (dashed 
black) . 

3.2.2. EoS predictions with the DBHF approach 

In Fig. |3l we show EoS predictions for symmetric matter (solid red) and neutron 
matter (dashed black) as obtained from the Idaho calculation described in the 
previous section. Equation (4) then gives the Eo S values for any a, a behavior 
which can be verified to be approximately true!^ 

The EoS from DBHF can be characterized as being moderately soft at low to 
medium density and fairly "stiff" at high densities. The predicted saturation density 
and energy for the symmetric matter EoS in Fig. [4] are equal to 0.185 fm -3 and - 
16.14 MeV, respectively, and the compression modulus is 252 MeV. For comparison, 
the same saturation obserbables as predicted by the other DBHF model presently 
on the markelP^are 0.181 fm -3 and -16.15 MeV for saturation density and energy, 
and 230 MeV for the incompressibility. 

The increased stiffness featured by the DBHF EoS at the higher densities 
originates from the strongly density-dependent repulsion inherent to the Dirac- 
Brueckner-Hartee-Fock method. In RefJ^, it is pointed out that constraints from 
neutron star phenomenology together with flow data from heavy-ion reactions sug- 
gest that such EoS behavior may be desirable. We will come back to this point 
later, in conjunction with neutron star predictions. 

At this point, it it useful to take a closer look at some of the predictions included 
in the analysis of such as relativistic mean field (RMF) models. Examples 

are those from RefsP^ES which use density-dependent ("DD") meson couplings 
and are fitted to the properties of nuclei up to about 0.15 fm -3 . They generate 
the steepest EoS and thus the largest pressures, see Fig. [5] An improvement to 
the traditional RMF description of nuclear matter can be obtained through the 
introduction of non- linear ( "NL" ) self-interactions of the a meson, such as done 
in the models of RefsfS3EH ; with the parametrization of 

RefP including both the 

S and the p mesons in the isovector channel. The corresponding EoS are much 
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less repulsive than those ol "DD" models (although the symmetry energy becomes 
very large at high density, possibly due to the absence of non-linearity and density 
dependence at the isovector level). 

Clearly, the pressure as a function of density plays the crucial role in building 
the structure of a neutron star. In Fig.[6]we show the pressure in symmetric matter 
as predicted by the Idaho calculation compared with constraints obtained from flow 
datafSSl The predictions are seen to fall just on the high side of the constraints and 
grow rather steep at high density. Comparing with Fig. [5j we see that the Idaho 
predictions are well below those of DD-RMF models at low to moderate density 
but nearly catch up with them at very high density, a description that would also 
be appropriate for the DBHF predictions of Ref!^ (red curve in Fig. [5}. Of all 
the cases studied in RefP^H DD-RMF models predict the largest maximum masses 
and radii and the lowest central densities. Thus, an equation of state where high 
pressure is sustained for a longer radial distance (moving away from the center 
of the star) will allow the maximum mass star to be heavier, larger, and more 
"diffuse" at the center. On the other hand, microscopic relativistic models, (such 

lonl 

as the DBHF calculation of Ref!^ or the present one, which are in fair agreement 
with each other), display a rather different density dependence of the pressure and 
produce smaller and more compact maximum mass stars. (All other EoS considered 
in RefP are softer and generate smaller maximum masses with smaller radii and 
larger central densities.) 

To conclude this section, we show in Fig. [7] the pressure in neutron matter (red 
curve) and /3-equilibrated matter (gr een) as predicted by DBHF calculations. The 
pressure contour is again from RefP^and was obtained from flow data together with 
the assumption of strong density dependence in the asymmetry term (indicated as 
"Asy_ stiff" in ReiP% 

3.2.3. The symmetry energy and related observables 

Here we will focus specifically on the symmetry energy and its impact on the struc- 
ture of neutron-rich nuclei, in particular the neutron skin thickness. 

In Fig. El we display the Idaho DBHF prediction for the symmetry energy by the 
solid red curve. The curve is seen to grow at a lesser rate with increasing density, an 
indication that, at large density, repulsion in the symmetric matter EoS increases 
more rapidly than in the neutron matter EoS. This can be understood in terms of 
increased repulsion in isospin zero partial waves (absent from neutron matter) as a 
function of density. Our predicted value for the symmetry pressure, L, (see Eq. ([38]) 
below), is 69.6 MeV. 

The various black dashed curves in Fig. |4] are obtained with the simple 
parametrization 

e sy m = C(p/po) J , (36) 
with 7 increasing from 0.7 to 1.0 in steps of 0.1, and C « 32 MeV. It seems that 
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Fig. 5. Pressure in symmetric matter predicetd by various models. The shaded area corre- 
sponds to the region of pressure consistent with the flow data analysed in Ref. 56 . [Figure 
reprinted with permission from T. Klahn. 51 Copyright (2006) by the American Physical Soci- 
ety. ^^JJ^^mg/^^^JP^Jy74/^J^35m^ 
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Fig. 6. Pressure in symmetric matter from the Idaho DBHF calculation. The shaded area corre- 
sponds to the region of pressure consistent with the flow data analysed in Ref. 56 . 



a value of 7 close to 0.8 gives a reasonable description of the DBHF predictions, 
although the use of different functions in different density regions would be best for 
an optimal ht. 
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I 1 , 1 , 1 , 1 
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P (fm 3 ) 

Fig. 7. Pressure in neutron (red curve) and baryon-lcpton (green curve) matter from the Idaho 
DBHF calculation. The shaded area corresponds to the region of pressure consistent with flow 
data and the inclusion of strong density dependence in the asymmetry terms. 56 



Considering that all of t he d ashed curves are commonly used parametrizations 
suggested by heavy-ion dataj^3 Fig. 2] clearly reflects our limited knowledge of the 
symmetry energy, particularly, but not exclusively, at the larger densities. 

As already mentioned in Sec. 2, from the experimental side intense effort is go- 
ing on to obtain reliable empirical information for the less known aspects of the 
EoS. Heavy- ion reactions are a popular way to seek constraints on the symmetry 
energy, through analyses of observables that are sensitive to the pressure gradient 
between nuclear and neutron matter. Isospin diffusion data from heavy-ion col- 
lisions, together with analyses based on isospin-dependent transport models, can 
provide information on the slope of the symmetry energy. 

Concerning the lower densities, isospin-sensitive observables can also be identi- 
fied among the properties of normal nuclei. The neutron skin of neutron-rich nuclei 
is a powerful isovector observable, being sensitive to the slope of the symmetry 
energy, which determines to which extent neutrons will tend to spread outwards 
to form the skin. It is the purpose of this section to systematically examine and 
discuss the symmetry energy properties in microscopic models and the correspond- 
ing neutron skin predictions. We will take the skin of 208 Pb as our representative 
isovector "observable" . 

Parity- violating electron scattering experiments are now a realistic option to de- 
termine neutron distributions with unprecedented accuracy. The neutron radius of 
208 Pb is expected to be measured within 0.05 fm thanks to the electroweak program 
at the Jefferson Laboratory!^ This level of accuracy could not be achieved with 
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Fig. 8. DBHF prediction for the symmetry energy (solid red) compared with various phenomeno- 
logical parametrizations (dashed black). See text for details. 



hadronic scattering. Parity- violating electron scattering at low momentum transfer 
is especially suitable to probe neutron densities, as the Z° boson couples primarily 
to neutrons. With the success of this program, reliable empirical information on 
neutron skins will be able to provide, in turn, much needed independent constraint 
on the density dependence of the symmetry energy. 

We form an energy functional based on the semi-empirical mass formula, where 
the volume and symmetry terms are contained in the isospin-asymmetric equation 
of state. Thus, we write the energy of a (spherical) nucleus as 

E(Z,A) = J d 3 r e(p(r),a(r))p(r) + 

+ J d 3 rf (\S7p(r)\ 2 + Pl {r)\ 2 ) + Coulomb term. (37) 

In the above equation, p and pi are the usual isoscalar and isovector densities, 
given by p n + p p and p n — p p , respectively, in terms of neutron and proton densities. 
a is the neutron asymmetry parameter, a — pi/p, and e(p, a) is the energy per 
particle in isospin-asymmetric nuclear matter. The latest Idaho Eo£p^will be used 
in Eq. (j37]), along with others from different models. 

From fits to nuclear binding energies, the constant /o in Eq. (|37|) is approxi- 
mately 70 MeV fm 5 , whereas the contribution of the term proportional to f3 was 
found to be minorPl Thus it is reasonable to neglect it. (The magnitude of (3 was 
estimated to be about 1 /4 in RefP, where it was observed that, even with variations 
of /? between -1 and +1, the effect of the j3 term on the neutron skin was negligibly 
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Fig. 9. (color online) Neutron (red) and proton (black) point densities as obtained from the DBHF 
model. 

small.) 

The parameters of the proton and neutron densities are obtained by minimizing 
the value of the energy functional, Eq. Q37p. assuming Thomas- Fermi distribution 
functions. Although simple, this method has the advantage of allowing a very direct 
connection between the EoS and the properties of finite nuclei. (It could be used, for 
instance, to determine a semi-phenomenological EoS by fitting it to both binding 
energies and charge radii of closed-shell nuclei.) 

In Fig. HI we show the proton and neutron Thomas- Fermi distributions for 208 Pb 
as obtained with the method described above and the DBHF model for the EoS. 
The predicted proton and neutron root-mean-square radii are 5.39 fm and 5.56 fm, 
respectively. 

Typically, predictions of the symmetry energy at saturation density encountered 
in the literature are in reasonable agreement with one another, ranging approxi- 
mately from 26 to 35 MeV. On the other hand, the slope of the symmetry energy 
is more model dependent and less constrained. This is seen through the symmetry 
pressure, defined as 

T o ( de sym{p) \ _ „ ( de n .m.{p) \ , Qa , 

Thus, L is sensitive to the gradient of the energy per particle in neutron matter 
(e n .m.)- As to be expected on physical grounds, the neutron skin, given by 

S = V< r 2 n > - ^j<rl > , (39) 

is highly sensitive to the same gradient. 

Values of L are reported to range from -50 to 100 MeV as seen, for instance, 
through the numerous parametrizations of Skyrme interactions (see RefP^ and ref- 
erences therein) , all chosen to fit the binding energies and the charge radii of a large 
number of nuclei. Heavy-ion data impose boundaries for L at 85 ± 25 

MeV,™ 
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with more stringent constraints being presently extracted!^ Other report state 
the constraints as 85 ±55 MeV. Also, a n early linear correlation is observed between 
the neutron skin S and the L parameter!^ 

Such phcnomcnological studies are extremely useful, but, ultimately, they must 
be compared with ab initio approaches in order to get true physical insight. The 
direct connection with the underlying nuclear forces will then facilitate the physical 
understanding, when combined with reliable constraints. 

As explained in Sec. 3.2.1, the DBHF model does not include three-body forces 
explicitcly, but effectively incorporates the class of TBF originating from the pres- 
ence of nucleons and antinucleons (the "Z-diagrams" in Fig. [2]), see previous dis- 
cussion. As the other main inpu t of our comparison, we will take the EoS from 
the microscopic approach of RefP. There (and in previous work by the same 
authors), the Brueckncr-Hartree-Fock (BHF) formalism is employed along with mi- 
croscopic three-body forces. In particular, in RefP^l the meson-exchange TBF arc 
constructed applying the same parameters as use d in the corresponding nucleon- 
nucleon (NN) potentials, which are: Argonne Vld§31 (V18), Bonn # (BOB), Ni- 
jmegen 93^2 (N93). The popular (but phenomenological) Urbana TBF^Sl (UIX) 
is also utilized in RefPS. Convenient parametrizations in terms of simple analytic 
functions are given in all cases and we will use those to generate the various EoS. We 
will refer to this approach, generally, as "BHF + TBF" . This comparison will bring 
up the discussion on microscopic models and explicit TBF outlined in Sec. 3.2. 

In Fig. 1101 we display Idaho DBHF predictions for the symmetry energy, solid 
black curve, along with those from V18, BOB, UIX, and N93. All values of the 
symmetry energy at the respective saturation densities are between 29 and 34 MeV. 
A larger spreading is seen in the L parameter, see left panel of Fig. [TTJ where the 
values range from about 70 to 106 MeV. The respective neutron skin predictions 
are shown on the vertical axis. 

We notice that all BHF+TBF models predict larger L, and thus larger neutron 
skins, as compared with DBHF, corresponding to a faster growth of the energy 
per particle in neutron matter relative to symmetric matter. This can be seen from 
Fig.rXUJ especially for the higher densities. The present calculations reveal that there 
are more subtle, but significant differences at low to medium densities as well. 

Some comments are in place concerning the nature of the L vs. S correlation. 
If a family of models is constructed so that they differ mostly in the slope of the 
symmetry energy, a linear correlation may be expected among L and the skinful 
The models we are considering, however, differ from one another in more than just 
the slope of the symmetry energy. Thus, although the general pattern is that a 
larger L corresponds to a larger skin, the relation is more complex, as the r.m.s. 
radius of the neutron distribution will, to some extent, receive feedback from the 
smaller or larger degree of attraction that binds neutrons and protons, and which 
depends on the symmetric matter EoS. 
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Fig. 10. (color online) Predictions for the symmetry energy from DBHF and various "BHF + 
TBF" models discussed in the text. 
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Fig. 11. Left panel: neutron skin of 208 Pb vs. the symmetry pressure for the models considered in 
the text; Right panel: neutron skin of 208 Pb vs. the curvature of the symmetry energy, K sym , for 
the same models. 



The next term in the expansion of the symmetry energy is the K sym parameter, 

which is a measure for the curvature of the symmetry energy. The neutron skin 
vs. K sym is shown in the right panel of Fig. [TT] for the various models. Although 
the values of K sym appear more spread out, the large negative values obtained 
with some of the parametrizations of the Skyrme model are not present. Those 
large negative values (as low as -600 MeV) produced by Skyrme models indicate 
a strongly downward curvature of the symmetry energy already at low to medium 
densities. We also notice from Fig. [TTJ (right panel) that no clear correlation can be 
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Fig. 12. Left panel: neutron skin of 208 Pb vs. the radius of a 1.4Mq neutron star for the models 
considered in the text. Right panel: star radius vs. the central energy density for the same models. 



identified between K sym and the neutron skin. 

We conclude this section by showing in Fig.[TS]the relation between the neutron 
skin of 208 Pb and the radius of a 1.4M Q neutron star for the models we are consid- 
ering. (A more detailed discussion of neutron stars and the EoS will be presented in 
Sec. 4.) Stellar matter contains neutrons in j3 equilibrium with protons, electrons, 
and muons. Tabulated values for the latest Idaho DBHF EoS can be found later in 
this a rtic le. We have applied /3-stability in the same way to all the various models 
of Refill starting from the given parametrized versions of the respective symmetric 
matter and neutron matter EoS. For the case of leptons, typically treated as a gas 
of non-interacting fermions, the equations of /3-stability (which amount to imposing 
energy conservation and charge neutrality), are elementary to solve (see Sec. 4.2). 

At subnuclear densities, all the EoS considered here are joined with the crustal 
equations of state from the work of Harrison and Wheeled (for energy densities 
between 10 and 10 11 g cm~ 3 ) and the work of Negele and Vautherir!§£l (for energy 
densities less than 1.7xl0 13 g cm -3 ). 

At first, Fig. 1121 can appear surprising, since it shows that larger skin does not 
necessarily imply larger radius. On the other hand, one must keep in mind that, as 
we argued previously, these models differ from one another in more than just the 
slope of the symmetry energy. Furthermore, the star radius probes higher densities 
than the skin does. Generally, if a family of models is constructed so that they differ 
mostly in the slope of the symmetry energy, a linear correlation is expected among 
L, the skin, and i?!^ 

In summary: Most models agree on the value of the symmetry energy around 
the saturation point within a few MeV, but we are far from a reasonable agreement 
on its derivative and even farther from agreement on its curvature. Although mi- 
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Fig. 13. The neutron (dashed line) and proton (dotted) single-particle potentials as a function of 
the momentum at fixed average density and a=0.4. The solid curve in the middle is the single- 
nucleon potential at o=0. 

croscopic models do not display as much spread as phenomenological ones, a point 
that comes out clearly from the present study is that a measurement of the neutron 
skin of 208 Pb with an accuracy of 0.05 fm, as it has been announced}^ should 
definitely be able to discriminate among the EoS from these microscopic models. 

3.2.4. Neutron and proton single-particle properties and the symmetry 
potential 

In this section we will concentrate specifically on predictions of isovector single- 
particle properties, such as the neutron and proton single-particle potentials and the 
closely related symmetry potential. These "observables" depend sensitively on the 
difference between neutron and proton properties in asymmetric matter and play 
an important role in simulations of heavy-ion collisions with neutron-rich nuclei. 

In terms of the G-matrix, calculated as described in Sec. 3.2.1, we write the 
single-nucleon potential (in the case of unequal Fermi levels for protons and neu- 
trons), as 

Ui(k) = Re[ < kq\G in \kq - qk > + ^ < kq\G lp \kq - qk >] , (41) 
<?<fcj q<k p F 

where i = n/p for neutron/proton, and k refers to states below and above the 
Fermi momentum (consistent with the "continuous choice" for the single-particle 
spectrum) . 

We begin by examining the momentum dependence of U n / P , the single neu- 
tron/proton potential in neutron-rich matter. In Fig. I13[ we show U n / p as a func- 
tion of the momentum and a fixed value of the asymmetry parameter, a = 0.4. The 
middle curve shows U n =U p at a =0. The total nucleon density considered in the 
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figure is equal to 0.185 fm -3 and corresponds to a Fermi momentum of 1.4 fm _1 , 
which is very close to our predicted saturation density. 

For increasing values of a (obviously, U n — U p for a—0), the proton potential 
becomes more attractive while the opposite tendency is observed in U n . This reflects 
the fact that the proton-neutron interaction, the one predominantly felt by the 
single proton as the proton density is depleted, is more attractive than the one 
between identical nucleons. Also, as it appears reasonable, the dependence on a 
becomes weaker at larger momenta. 

The role of the momentum dependence of the symmetry potential in heavy- 
ion collisions was examined and found to be important. Symmetry potentials with 
and without momentum dependence can lead to significantly different predictions 
of collision observa blesj^ More recently, these issues have been revisited in hot 
asymmetric matter 

Regarding U n / p as functions of the asymmetry parameter a, one can easily verify 
that the following approximate relation applies: 

U n/p {k, k F ,a) w U n / P (k, k F , a = 0) ± U sym {k, k F )a , (42) 

with the ± referring to neutron/proton, respectively. Figure [14] displays the left- 
hand side of Eq. (|42|) for fixed density and nucleon momentum and clearly reveals 
the linear behaviour of U n / p as a function of a. Thus, one can expect isospin splitting 
of the single-particle potentials to be effective in separating the collision dynamics 
of neutrons and protons. Although the main focus of Fig. [14] is the a dependence, 
predictions are displayed for three different potentials, Bonn Ap^Bp^and cPThese 
three models differ mainly in the strength of the tensor force, which is mostly carried 
by partial waves with isospin equal to and thus should fade away in the single- 
neutron potential as the neutron fraction increases. In fact, the figure demonstrates 
reduced differences among the values of U n predicted with the three potentials at 
large a. 

Already several decades ago, it was pointed out that the real part of the nuclear 
optical potential depends on the asymmetry parameter as in Eq. (142 j) !^ Then, the 
quantity 



Un + Up 



= U , (43) 



which is obviously the single-nucleon potential in absence of asymmetry, should 
be a reasonable approximation to the isoscalar part of the optical potential in the 
interior of a nucleus. The momentum dependence of Uq (which is shown in Fig. [13] 
as the a=0 curve), is important for extrac ting information about the symmetric 
matter EoS and is reasonably agreed upon I 75 | 76 | 77 | 78 | 79 | 80 | 81 | 82 | 
On the other hand, 



En Hp - Tj (U) 

2 a ~ u sym 
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Fig. 14. The neutron and proton single-particle potentials as a function of the asymmetry pa- 
rameter at fixed average density and momentum equal to the average Fermi momentum. 



should be comparable with the Lane potential, or the isovector part of the nuclear 
optical potential. (In the two equations above the dependence upon density, mo- 
mentum, and asymmetry has been suppressed for simplicity.) We have calculated 
U sym as a function of the momentum, or rather the corresponding kinetic energy. 
The predictions obtained with Bonn B are shown in Fig. 1151 They are compared 
with the phenomenological expression^] 

U Lane = a-bT, (45) 

where T is the kinetic energy, a » 22 - MMeV, b « 0.1 - 0.2MeV. 

We observe that the strength of the (DBHF) predicted symmetry potential 
decreases with energy, a behavior w hich is consistent with the empirical information. 
The same comparison is done in Ref 1^ starting from a phenomenological formalism 
for the single-nucleon potential! ^^ There, it is shown that it is possible to choose 
two sets of parameters which lead to similar values of the symmetry energy but 
exactly opposite tendencies in the energy dependence of the symmetry potential 
as well as opposite sign of the proton-neutron mass splitting. As a consequence of 
that, these two sets of parameters lead to very differentpredictions for observables 

lo r I 

in heavy-ion collisions induced by neutron-rich nuclei."- This fact suggests that 
constraints from "differential" observables, namely those specifically sensitive to 
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Fig. 15. The symmetry potential close to saturation density. 




Fig. 16. The neutron and proton effective mass as a function of the asymmetry parameter at fixed 
average density. 



the difference between proton and neutron properties in asymmetric matter, are 
very much needed. 

The effective masses for proton and neutron corresponding to the single-nucleon 
potentials of Fig. [T3] are shown in Fig.[T()]as a function of a. The predicted effective 
mass of the neutron being larger tha n the proton's is a trend shared with micro- 
scopic non-relativistic calculations !^ In the non-relativistic case, one can show from 
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Table 1. np total effective cross sections in symmetric matter calculated 
with various many-body models, see text for details. 



kpUm- 1 ) qo(McV) o° p BHF {mb) o B p HF+TBF (mb) u BBF (rab) 



1.1 
1.1 
1.1 
1.4 
1.4 
1.4 
1.7 
1.7 
1.7 



250 
300 
350 
250 
300 
350 
250 
300 
350 



34.38 
23.14 
20.63 
26.74 
17.26 
16.77 
17.20 
15.06 
12.33 



44.65 
29.01 
23.56 
31.25 
25.28 
21.17 
19.03 
17.59 
13.99 



51.74 
31.85 
25.62 
39.82 
30.96 
25.23 
29.12 
25.02 
21.04 



elementary arguments based on the curvature of the (parabolic) single-particle po- 
tential that a more attractive potential, as the one of the proton, leads to a smaller 
effective mass. In the DBHF quasi-particle approximation described in Sec. 3.2.1, 
one assumes momentum-independent nucleon self-energies, Us and Uy, with a van- 
ishing spacial component of the vector part. In such limit, the "Dirac mass" and 
the non-relativistic or "Landau mass" should coi ncid e to leading order in the ex- 
pansion of the relativistic single-particle potential!*^ It is therefore not surprising 
that our proton and neutron effective masses, although defined as Dirac masses 
(m* = m + Us,i), display a trend that is similar to the one generally expected in 
the non-relativistic case. Most likely, the main reason for the discrepancy in the pre- 
dicted isospin mass splitting between the Idaho DBHF calculation and the DBHF 
model from RefP 

is in the way the scalar and vector fields are treated, particularly 
their momentum dependence. 

Closely related to the effective mass is the concept of in-medium effective cross 
section. Transport equations, such as the Boltzmann-Uehling-Uhlenbeck equation, 
describe the evolution of a non-equilibrium gas of strongly interacting hadrons 
drifting in the presence of the mean field while undergoing two-body collisions. 
Thus heavy ion collision simulations require the knowledge of in-medium two-body 
cross sections as well as the mean field. Both should be calculated microscopically 
and self-consistently starting from the basic nuclear forced^ and are closely related 
to the EoS. 

To be fully applicable throughout the evolution of the colliding system, where 
high densities may be involved, in-medium cross sections must go beyond the con- 
ventional BHF. Although an extensive discussion of in- medium cross sections will 
not be presented here, we show in Tables 1 and 2 a comparison of the latest Idaho 
(DBHF) in-medium cross sections and those obtained from the BHF+TBF model 
of Lombardo et aZ.P^l 

We recall that the relativistic elastic differential cross section is calculated from 



a{6) 




\G(e)\\ 



(46) 



ATT 2 fl 4 S* 
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Table 2. As in the previous Table, for identical particles. 



kpUm' 1 ) 


qo(MeV) 




BHF+TBF i ,s 
a NN ( mb ) 




1.1 


250 


18.00 


18.15 


22.98 


1.1 


300 


16.41 


17.47 


19.48 


1.1 


350 


17.08 


18.55 


19.51 


1.4 


250 


15.72 


13.74 


19.76 


1.4 


300 


13.70 


14.43 


17.23 


1.4 


350 


16.31 


16.89 


19.46 


1.7 


250 


18.05 


7.87 


13.36 


1.7 


300 


17.93 


9.98 


13.98 


1.7 


350 


13.96 


11.48 


16.43 



where s* = 4((to*) 2 + q 2 ). Notice that the relativistic amplitude, G, is related to 
the non-relativistic one through G = |p-G„. r ., with E* = ((m*) 2 + g 2 ) 1 / 2 . 

In Tables 1 and 2, qo signifies the momentum of either nucleon in their center-of- 
mass frame, for simplicity assumed to coincide with the nuclear matter rest frame. 
The 4th and 5th columns show the results obtained by Lombardo et al. with BHF 
calculations implemented or not with TBF. There is generally a fair amount of 
agreement between the DBHF and the BHF+TBF results. For the np case, energy 
and density dependence appear consistent among the two sets of results displayed 
in columns 3 and 4, although the BHF+TBF results are generally larger. For the 
case of identical nucleons, the DBHF values are in good agreement with those from 
BHF+TBF, with the exception of the highest densities. This is to be expected 
in view of the larger differences existing between the DBHF and the BHF+TBF 
predictions of the EoS at high density. 

Besides being a crucial part of the input for transport models, in-medium ef- 
fective cross sections are important in their own right as they allow to establish an 
immediate connection with the nucleon mean free path, A, one of the most funda- 
mental properties characterizing the propagation of nucleons through matter!^ It 
can be written as 

A = j , (47) 

&ppPp t &pnPn 

for a proton propagating through matter with proton and neutron densities equal 
to p p and p n , respectively. The mean free path enters the calculation of the nuclear 
transparency function, which is the probability that the projectile will pass through 
the target without interacting. This is closely related to the total reaction cross 
section of a nucleus, which can be used to extract nuclear r.m.s. radii within Glauber 
models!^ Therefore, microscopic in-medium isospin-dependent NN cross sections 
can ultimately help obtain information about the properties of exotic, neutron- 
rich nuclei, such as, for instance, the target density. The latter is a very useful 
information, especially for exotic nuclei or nuclei with large neutron skins (where 
densities cannot be probed with conventional charged lepton scattering). 
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4. Neutron Star Properties 

Proceeding in our review of the various systems and phenomena where the EoS 
plays a chief role, we will now take a look at some of the most exotic systems in 
the universe, compact stars. 

4.1. A brief review of neutron star structure and available 
constraints 

Neutron stars are stable configurations containing the most dense form of matter 
found in the universe. They are therefore unique laboratories to study the properties 
of highly compressed (cold) matter. Furthermore, the possibility of exploring the 
structure of neutron stars via gravitational waves^ makes these exotic objects even 
more exciting. 

The densities found in neutron stars range from the density of iron to several 
times normal nuclear density. In the low-density region, matter is highly compressed 
and fully ionized and consists of electrons and ions of iron. As density increases, 
charge neutrality requires matter to become more neutron rich. In this density range 
(10 7 < p < 10 11 g cm -3 ), neutron-rich nuclei appear, such as isotopes of 62 Ni, 86 Kr, 
84 Se, 80 Zn, 124 Mo, 122 Zr. Above densities of approximately 10 11 g cm -3 , free neu- 
trons begin to form a continuum of states. The inner crust is then a compressed 
solid with a fluid of neutrons. At densities equal to approximately 1/2 of saturation 
density, clusters begin to merge into a continuum. In this phase, matter is a uni- 
form fluid of neutrons, protons, and leptons. Protons and neutrons can conde nse 
in superfluids and superconducting states below temperatures of about 10 9 K.^U 
Above a few times nuclear matter density, the actual composition of stellar matter 
is not known. Strange baryons can appear when the nucleon chemical potential is 
of the order of their rest mass. We will address the issue of strange matter later 
in this article. At even higher densities, transitions to other phases are speculated, 
such as a deconfined, rather than hadronic, phase. The critical density for such 
transition cannot be predicted reliably because it lies in a range where QCD in non 
perturbative 

Recently, a semi- algebraic expression for the EoS of cold quark matter has been 
derived within the Dyson- Schwinger formalism!^ The Dyson-Schwinger framework 
can address both confinement and dynamical chiral symmetry breaking]^ unlike 
other models commonly used in conjunction with cold quark matter, such as the 
Nambu-Jona-Lasinio model or bag-type models. It will be interesting to see ap- 
plications of the model described in RefP 

to dense astrophysical systems. For 
that purpose, the quark matter EoS will need to be supplemented by the nuclear 
matter EoS. In the meantime, it is encouraging that the model predicts coincident 
deconfinement and chiral symmetry restoring transitions. 

The possibility has been speculated that the most stable state at zero pressure 
may be it, d, s quark matter instead of iron. This would imply that strange quark 
matter is the most stable (in fact, the absolutely stable) state of strongly interacting 
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matter, as originally proposed by Bodmer WittenpSl and Terazawa!^ In such 
case, hyperonic and hybrid stars would have to be metastable with respect to stars 
composed of stable three-flavor strange quark matter j^l which is lower in energy 
than two-flavor quark matter due to the extra Fermi levels open to strange quarks. 
Whether or not strange quark stars can give rise to pulsar glitches, (which are 
observed sudden small changes in the rotational frequency of a pulsar) , may be a 
decisive test of the strange quark matter hypothesis.^ 

The maximum mass and the radius of a neutron star are sensitive to different 
aspects of the EoS. The maximum mass is mostly determined by the stiffness of the 
EoS at densities greater than a few times saturation density. Non-nucleonic degrees 
of freedom, which typically make their appearance at those densities, are known to 
have a considerable impact on the maximum mass of the star. Causality imposes 
a limit of about 3 solar masses to the maximum mass, the existence of which 
is inherent to general relativity and is predicted by the equation of hydrostatic 
equilibrium, 

dP(r) G (P(r) + e(r))(M(r) + 47rr 3 P(r)/c 2 ) 



dr c 2 r{r -2GM{r)/c 2 ) 

with 

dM(r 



(48) 



4tt V) , (49) 
dr 

where e is the total mass-energy density. The pressure is related to the en- 
ergy /particle through 

^)=P 2 ^- (50) 

It's worth recalling that no mass limit exists in Newtonian gravitation. 

The star radius is mainly sensitive to the slope of the symmetry energy. In par- 
ticular, it is closely connected to the internal pressure (that is, the energy gradient) 
of matter at densities between about 1.5,00 and 2-3po-^ 

Lattimer and Prakash investigated the maximum central density question!^ 
Combining the causality limit (R >2.87GM/c 2 ) with the central density-mass re- 
lation implied by the Tolman VII solutionp 21 

p = p c (l - r 2 /R 2 ) , (51) 

they obtained the relation 

p c M 2 < 15.3 x 10 15 M| g cm- 3 . (52) 

The maximum gravitational mass and the corresponding radius are the typical 
observables used to constraint the EoS. Much more stringent constraints could be 
imposed on EoS if these two quantities could be determined independently from 
each other. At this time, the hea vies t neutron star (with accurately known mass) 
has a mass of 1.671 ± 0.008 M Q ,-^ but neutron stars with masses above 2 solar 
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masses are expected to exist JEH3 Mass limits are obtained from observations of 
binary systems, either a two-neutron star system or a neutron star and a ma ssive 
companion, such as a white dwarf. The pulsar in the Hulse- Taylor binary 

system™ 

has a mass of 1.4408± 0.0003 M©, to date the best mass determination. 

Measurements of the radius are considerably less precise than mass 
measurements.^3 The radiation or photospheric radius, i?oo, is related to the actual 
stellar radius by B, x = R{l-R/R s )~ 1/2 , where R s = 2GM/R is the Schwarzschild 
radius. Estimates are usually based on thermal emission of cooling stars, including 
rcdshifts, and the properties of sources with bursts or thermonuclear explosions at 
the surface. 

Neutrinos are also used as a probe of the EoS. Neutrino luminosity is controlled 
by several factors including the total mass of the (proton-neutron) star and the 
opacity of neutrinos at high densities, which of course is sensitive to the EoS of 
dense matter. 

Gravitational wave 

s 90 

are a less conventional way to probe neutron star proper- 
ties. Emission of gravitational waves causes orbital decay, the observation of which 
would allow an estimate of the moment of inertia and thus, together with an ac- 
curate measurement of the mass, impose stronger boundaries on the EoS. A mea- 
surement of the moment of inertia within 10%, together with the information on 
the mass, would be able to discriminate among various 

EosEZl To date the best 



The minumum mass of a neutron star is also a par ame ter of interest. For a 

1071 

cold, stable system, the minimum mass is about 0.09 MQ^-'The smallest reliably 
estimated neutron star mass is the com panion of the binary pulsar J 1756-2251, 
which has a mass of 1.18± 0.02 M©PSI 

4.2. Comparing predictions of microscopic models 

As constraints from nuclear physics and/or astrophysics promise to become more 
stringent, it is important to understand and compare how the nature of the various 
predictions is related to the features of each model. In microscopic approaches, the 
tight connection with the underlying forces facilitates physical interpretation of the 
predictions in terms of the characteristics of the nuclear force and its behavior in 
the medium. Motivated by these considerations, in this section we calculate several 
neutron star properties, for static and/or rotating stars, using equations of state 
based on different microscopic models. These will be the same as those used in 
Sec. 3.2.3 to examine symmetry energy and neutron skin predictions, namely BOB, 
V18, N93, and UIX, along with Idaho DBHF. 

In Sec. 3.2.3, we compared the predictions of the neutron skin in 208 Pb by these 
models and correlated them with differences in the slope of the symmetry energy. 
Model differences become larger at high-density and will naturally impact neutron 
star predictions. 



determination of the moment of i 
would rule out only very soft EoS. 
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Table 3. DBHF equation of state of pure neutron matter. 



Baryon density(l/cm 3 ) Energy density (j/cm 3 ) Pressure(dj/ne/cm 2 ) 






145902E+38 





245257E+14 





434992E+32 





231688E+38 





389881E+14 





874118E+32 





345843E+38 





582634E+14 





156362E+33 





492421E+38 





830558E+14 





257664E+33 





675475E+38 





114075E+15 





468635E+33 





781946E+38 





132151E+15 





621198E+33 





899057E+38 





152061E+15 


o 


816935E+33 





102731E+39 





173905E+15 


o 


111487E+34 





116722E+39 





197791E+15 





156721E+34 





131929E+39 





223833E+15 





223210E+34 





148402E+39 





252158E+15 





319807E+34 





166192E+39 





282905E+15 





463315E+34 





185350E+39 





316247E+15 





681545E+34 





205927E+39 





352381E+15 





993009E+34 





227973E+39 





391534E+15 





141505E+35 





251538E+39 





433969E+15 





203050E+35 





276674E+39 





480078E+15 





293314E+35 





303432E+39 





530247E+15 





409031E+35 





331861E+39 





584918E+15 





558832E+35 





362012E+39 





644589E+15 


o 


746722E+35 





393937E+39 





709762E+15 





965629E+35 





427685E+39 





780907E+15 





122694E+36 





463308E+39 





858621E+15 





152765E+36 


o 


500856E+39 





943372E+15 


o 


185795E+36 





540380E+39 





103570E+16 





223317E+36 





581930E+39 





113630E+16 





268106E+36 





625557E+39 





124615E+16 





324163E+36 





671312E+39 





136650E+16 





395072E+36 





719245E+39 





149890E+16 





485417E+36 





769408E+39 





164502E+16 





589746E+36 





821850E+39 





180610E+16 





707244E+36 





876622E+39 





198388E+16 





846053E+36 





933776E+39 





218030E+16 





100924E+37 





993361E+39 





239765E+16 





120127E+37 





105543E+40 





263842E+16 





142307E+37 





112003E+40 





290521E+16 





167838E+37 





118721E+40 





320102E+16 





197388E+37 





125703E+40 





352921E+16 





231442E+37 





132954E+40 





389345E+16 





270588E+37 





140478E+40 





429785E+16 





315475E+37 





148280E+40 





474691E+16 





366826E+37 





156366E+40 





524560E+16 





425495E+37 





164741E+40 





579940E+16 





492059E+37 





173410E+40 





641418E+16 





567978E+37 





182378E+40 





709679E+16 





654219E+37 



As mentioned in Sec. 3.2.3, at subnuclear densities all the EoS considered here 
are joined with the crustal equations of state from Harrison and WheeleJ^ and 
Negele and VautherinP^ The composition of the crust is crystalline, with ligh 
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Table 4. As in the previous Table, for /3-equilibrated matter. 



Baryon density(l/cm 3 ) Energy density (j/cm 3 ) Pressure(dj/ne/cm 2 ) 






145902E+38 





245236E+14 





416322E+32 





231688E+38 





389827E+14 





822804E+32 





345843E+38 





582501E+14 





142104E+33 





492421E+38 





830269E+14 





230061E+33 





675475E+38 





114019E+15 





409814E+33 





781946E+38 





132073E+15 





540126E+33 





899057E+38 





151955E+15 





711724E+33 





102731E+39 





173765E+15 


o 


970498E+33 





116722E+39 





197605E+15 





137046E+34 





131929E+39 





223589E+15 





201627E+34 





148402E+39 





251837E+15 


o 


261928E+34 





166192E+39 





282481E+15 





330297E+34 





185350E+39 





315683E+15 





517205E+34 





205927E+39 





351619E+15 





757028E+34 





227973E+39 





390495E+15 


o 


110037E+35 





251538E+39 





432555E+15 





160843E+35 





276674E+39 





478138E+15 





235092E+35 





303432E+39 





527601E+15 





334557E+35 





331861E+39 





581355E+15 


o 


465834E+35 





362012E+39 





639861E+15 





632143E+35 





393937E+39 





703594E+15 





832228E+35 





427685E+39 





773030E+15 





107139E+36 





463308E+39 





848703E+15 





134775E+36 


o 


500856E+39 





931088E+15 





165606E+36 





540380E+39 





102069E+16 





200271E+36 





581930E+39 





111814E+16 





241417E+36 





625557E+39 





122434E+16 





293125E+36 





671312E+39 





134047E+16 





358691E+36 





719245E+39 





146797E+16 





442173E+36 





769408E+39 





160841E+16 





539098E+36 





821850E+39 





176298E+16 





649283E+36 





876622E+39 





193337E+16 





779667E+36 





933776E+39 





212141E+16 





932987E+36 





993361E+39 





232930E+16 





111395E+37 





105543E+40 





255941E+16 





132452E+37 





112003E+40 





281437E+16 





156916E+37 





118721E+40 





309721E+16 





185697E+37 





125703E+40 





341144E+16 





219046E+37 





132954E+40 





376060E+16 





257227E+37 





140478E+40 





414881E+16 





301377E+37 





148280E+40 





458067E+16 





352202E+37 





156366E+40 





506130E+16 





410655E+37 





164741E+40 





559636E+16 





477543E+37 





173410E+40 





619202E+16 





554149E+37 





182378E+40 





685522E+16 





641233E+37 



or heavy^ metals and electron gas. The DBHF equations of state for neutron 
matter and for /3-equilibrated matter are given in Table 3 and Table 4, respectively. 
Some points of the EoS used for the crust are displayed in Table 5. 
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Table 5. Equation of state used for the crust. See text for details. 



Baryon dcnsity(l/cm 3 ) Energy density (g/cm?) Pressure(dj/ne/cm 2 ) 



o 


59701000E+25 





99998600E+01 





40721410E+12 


o 


19099900E+26 





31992190E+02 





47174340E+14 





38376300E+26 





64279760E+02 





33550870E+15 





15492600E+27 





25949850E+03 





88776630E+16 





31128200E+27 


o 


52139300E+03 





38114990E+17 


o 


62543700E+27 





10476010E+04 





15215920E+18 





25249000E+28 





42291870E+04 





20992040E+19 





50731100E+28 





84974150E+04 





74261220E+19 





20480200E+29 





34304300E+05 





87378890E+20 


o 


41149300E+29 


o 


68925200E+05 


o 


29327530E+21 





13164600E+30 





22050990E+06 





21592540E+22 





26450500E+30 


o 


44305570E+06 





70850350E+22 





53144400E+30 





89020260E+06 





23123900E+23 





21453100E+31 





35937570E+07 





19523170E+24 


o 


43102100E+31 





72207080E+07 


o 


54232090E+24 


o 


13787700E+32 





23100980E+08 





28166750E+25 





34951600E+32 


o 


58569360E+08 





10057550E+26 


o 


70215700E+32 


o 


11767960E+09 





25560490E+26 


o 


28334500E+33 


o 


47507410E+09 





15831900E+27 





56915000E+33 





95453350E+09 





38815460E+27 





22959600E+34 





38534730E+10 





22901680E+28 





46107400E+34 





77425110E+10 





55287600E+28 





14734400E+35 





24770450E+11 





23888140E+29 


o 


29580200E+35 


o 


49769610E+11 





57354730E+29 





59374600E+35 


o 


99998600E+11 


o 


13755860E+30 


o 


30685800E+36 





51813430E+12 


o 


66738040E+30 





60373600E+36 


o 


10203670E+13 


o 


12773240E+31 





10985200E+37 


o 


18580330E+13 


o 


20087290E+31 





15931600E+37 





26960450E+13 





26615840E+31 





19888300E+37 





33666280E+13 





32106350E+31 





23844100E+37 





40373370E+13 





38075420E+31 





26810600E+37 





45404400E+13 





42918160E+31 





29776600E+37 





50436140E+13 





48089980E+31 





32742200E+37 





55468420E+13 





53594420E+31 





35707300E+37 





60501410E+13 





59429390E+31 





38672000E+37 





65534760E+13 





65589600E+31 





41636500E+37 





70568810E+13 





72068480E+31 





44600500E+37 





75603230E+13 





78857870E+31 





48552100E+37 





82316740E+13 





88379610E+31 





51515000E+37 





87352400E+13 





95862580E+31 





54477600E+37 





92388600E+13 





10362880E+32 





57439900E+37 





97425150E+13 





11167080E+32 





60401800E+37 





10246220E+14 





11998160E+32 





65337700E+37 





11085860E+14 





13440860E+32 





69285800E+37 





11757640E+14 





14645200E+32 





73233400E+37 





12429500E+14 





15892440E+32 





77180500E+37 





13101440E+14 





17181270E+32 





82113400E+37 





13941490E+14 





18848650E+32 





86059200E+37 





14613590E+14 





20226200E+32 





89018100E+37 





15117700E+14 





21284120E+32 





93948300E+37 





15958020E+14 





23093300E+32 





98877600E+37 





16798420E+14 





24958070E+32 
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Fig. 17. (color online) Proton (solid black), electron (dashed red), and muon (dotted blue) fractions 
in /3-stable matter as a function of total baryon density as predicted by the DBHF model. 



The proton fraction in /3-stable matter is calculated by imposing energy con- 
servation and charge neutrality. The resulting algebraic equations can be found in 



1 0^ 

standard literature! mjoi -j^g contribution to the energy density from the electrons 
is written as 



he 

4^ 



whereas for muons we write 



p^m^c 2 + (Hc)" 



107r 2 m M c 2 



(53) 



(54) 



These contributions are added to the baryonic part to give the total energy den- 
sity. The derivative of the total energy/particle with respect to the fraction of a 
particular particle species is the chemical potential of that species. The conditions 



Pp + Me = Mn ; Mm = Me ] Pp = Pn + Pe 



(55) 



allow to solve for the densities (or fraction) of protons, electrons, and muons. Near 
the saturation density, when the muon f ract ion is close to zero, one can estimate 
the equilibrium proton fraction, x p , to 

beP 



'4e 



p ~ \ he J /{3n P ° ] 



(56) 



The fractions of protons, electrons, and muons as predicted by DBHF are shown 
in Fig. [TT] The critical density for the proton fraction to exceed approximately 1/9 
and, thus, allow cooling through the direct Urea processes, 



n — >■ p + e + v e and p - 
is about 0.36 - 0.39 far 3 . 



(57) 
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In Fig. 1181 we show the mass-radius relation for a sequence of static neutron 
stars as predicted by the various models considered previously, see Fig. [TJ3 All 
models besides DBHF share the same many-body approach (BHF+TBF) but differ 
in the two-body potential and TBF employed. The resulting differences can be much 
larger than those originating from the use of different many-body approaches. This 
can be seen by comparing the DBHF and BOB curves, both employing the Bonn 
B interaction (although in the latter case the non-reltivistic, r-space version of the 
potential is adopted). Overall, the maximum masses range from 1.8M@ (UIX) to 
2.5M (BOB). Radii are less sensitive to the EoS and range between 10 and 12 km 
for all models under consideration, DBHF or BHF+TBF. Concerning consistency 
with present constraints, an initial observation of a neutron-star-white dw arf binary 
system suggested a neutron star mass (PSR J0751+1807) of 2T±O.2M P22 Such 
observation would imply a considerable constraint on the high-density behavior 
of the EoS. On the other hand, a dramatically reduced value of 1.3±O.2M was 
reportedly later, which does not invalidate any of the theoretical models under 
consideration. 

The model dependence is shown in Fig. [19] for the case of rapidly rotating 
stars. The 716 H z frequency corresponds to the most rapidly rotating pulsar, PSR 
J1748-2446p23 although recently an X-ray burst oscillation at a frequency of 1122 
Hz has been reporte d 1 which may be due to the spin rate of a neutron star. 
As expected, the maximum mass and the (equatorial) radius become larger with 
increasing rotational frequency. 

Another bulk property of neutron stars is the moment of inertia, I. In Fig. |2"01 
we show the moment of inertia at different rotational speeds (again, for all models), 
whereas in Fig.[5T]we display the moment of inertia corresponding to the maximum 
mass at different rotational frequencies. These values are not in contradiction with 
observations of the Crab nebula luminosity. From that, a lower bound on the mo- 
ment of inertia was inferred to be / >4-8 x 10 44 g cm 2 , see 

RefP and references 

therein. The size of lM max changes from model to model in line with the size of the 
maximum mass, see Fig. [2T] 

Lastly, we calculate the gravitational redshift predicted by each model. The 
redshift is given by 



This simple formula can be derived considering a photon emitted at th e su rface of a 
neutron star and moving towards a detector located at large distance!^! The pho- 
ton frequency at the emitter (receiver) is the inverse of the proper time between two 
wave crests in the frame of the emitter (receiver). Assuming a static gravitational 
field, and writing goo as the metric tensor component at the surface of a nonrotat- 
ing star yield the equation above. Naturally the rotation of the star modifies the 

a All neutron star properties are calculated from public software downloaded from the website 
http://www.gravity.phys.uwm.edu/rns 




(58) 
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R (km) 

Fig. 18. (color online) Static neutron star mass-radius relation for the models considered in the 
text. 

metric, and in that case different considerations need to be applied which result in 
a frequency dependence of the redshift. We will not consider the general case here. 

From Fig. [521 it appears that the gravitational redshift is not very EoS- 
dependent (compare, for instance, the values at the maximum mass for each model), 
an indication that the gravitational profile at the surface of the star is similar in all 
models. Thus measurements of z may not be the best way to discriminate among 
different EoS. 

Clearly, at the densities probed by the interior of neutron stars the model de- 
pendence is large, but presently available constraints are still insufficient to dis- 
criminate among these EoS. At very high density (p between five and ten times 
normal density) , the most repulsive symmetric matter energies are produced with 
BOB, V18, DBHF, N93 and UIX, in that order. This explains the maximum mass 
predictions, which depend mostly on the absolute repulsion present in the symmet- 
ric matter EoS. In pure neutron matter, N93 follows right after BOB (again, from 
largest to smallest repulsion). This indicates a somewhat different balance of attrac- 
tion/repulsion when only T=l contributions are included. The symmetry energy, 
see Fig. I10| which depends entirely on the repulsion of neutron matter relative to 
symmetric matter and whose density dependence controls observables such as the 
neutron skin, is largest in N93, followed by BOB, V18, UIX, and DBHF. 

The model dependence we observe comes from two sources, the two-body po- 
tential and the many-body approach, specifically the presence of explicit TBF or 
Dirac effects. The dependence on the two-body potential is very large. Typically, 
the main source of model dependence among NN potentials is found in the strength 
of the tensor force. Of course, differences at the two-body level impact the TBF as 
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Fig. 19. (color online) Mass-radius relation for the models considered in the text and for different 
rotational frequencies. 



well, whether they are microscopic or phenomenological (as in the case of UIX). 

On the other hand, when comparing DBHF and BOB, we are looking at differ- 
ences stemming from the many-body scheme, as the two models share the same NN 
potential. In the BOB model, repulsion grows at a much faster rate than in DBHF, 
and more strongly so in neutron matter. (Hence, the much larger symmetry energy 
with BOB). As an example, at about 6 times normal density the DBHF energy of 
symmetric matter is 67% of the BOB energy and only 51% at ten times normal 
densities. In pure neutron matter, those ratios become 49% and 30%, respectively. 
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Fig. 20. (color online) Moment of inertia for the models considered in the text and for different 
rotational frequencies. 

Thus, the inclusion of the microscopic TBF in the BHF model introduces consider- 
able more repulsion than the Dirac effects through highly non-linear termsPSl Both 
attractive and repulsive TBF are required for a realistic description of the satura- 
tion point. The density dependence of the repulsive terms is obviously stronger and 
thus dominates at high density. Furthermore, it appears that this is especially true 
in neutron matter. 

In Ref.ESI it is shown that the largest contribution to the net TBF originates 
from Z-diagrams such as shown in Fig. 2, a fact which confirms the validity of the 



January 31, 2010 19:6 WSPC/INSTRUCTION FILE Sammarruca'ijmpe 



The Microscopic Approach to Nuclear Matter and Neutron Star Matter 39 



4.0 

E 
o 

□> 

% 3.2 

eg 
t 

CD 

£ 2.4 

M — 

O 



E 1.6 

o 



■ DBHF 








• BOB 








- A N93 










♦ UIX 








• 


★ V18 




• 


• 




• 


• 
















* 


★ 


★ 


* 


t 




■ 


■ 






▲ 


▲ 


▲ 


▲ 


▲ 












♦ 


♦ 


♦ 


♦ 


♦ 





300 



1200 



600 900 
Frequency (Hz) 

Fig. 21. (color online) Moment of inertia corresponding to the maximum mass for the models 
considered in the text as a function of the rotational frequency. 
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Fig. 22. (color online) Gravitational redshift for all models. For each model, the corresponding 
sequence of static stars is considered. 



DBHF approximation. The actual amount of repulsion, however, seems to depend 
sensitively upon the TBF parametrization. One of the advantages of the DBHF 
method is the ability to avoid possible inconsistencies between the parameters of 
the two- and the three-body systems. 

At this time, available constraints cannot pin down the high-density behavior of 
the EoS. Nevertheless, we argue again that microscopic models allow for a deeper 
insight into the origin of the observed physical effects, and should be pursued along 
with improved constraints. 



January 31, 2010 19:6 WSPC/INSTRUCTION FILE Sammarruca'ijmpe 



40 Francesco, Sammarruca 

5. Non-Nucleonic Degrees of Freedom 
5.1. Introduction 

At densities close to normal nuclear density, protons and neutrons are the only 
baryonic degrees of freedom. As density increases, other baryons begin to appear, 
such as strange baryons or isospin 3/2 nucleon resonances. Hyperonic states can be 
classified according to the irreducible representation of the SU(3) group. The octet 
of baryons that can appear in neutron matter includes nucleons, A, £ 0,± , and S 0,_ . 

Neglecting the nucleon-hyperon interaction, the threshold for stable hyperons 
to exist in matter is determined by comparing the hyperon mass with the neutron 
Fermi energy, which is the largest available energy scale in neutron-rich matter. We 
consider cold neutron stars, after neutrinos have escaped. Strange baryons appear 
at about 2-3 times normal density,^ an estimate which is essentially model inde- 
pendent, through the processes n + n — > and n + n — > n + A. The equilibrium 
conditions for these reactions are 



2/i„ = Mp + Ms- 



Also, we have 

pe — Pu j Pn 

the equations above being special cases of 



; Pn = Pa ■ (59) 
= fi p + p e , (60) 
- IPe , (61) 



where b and q are the baryon number and the charge (in units of the electron 
charge) of the particular species with chemical potential p. Together with the charge 
neutrality condition and baryon number conservation, 

Pp = Pe+ Pu + PS- ; P = Pn + Pp + PE- + PA , (62) 

the above system allows to determine the various particle fractions. 

Clearly, the composition of matter at supra-nuclear densities determines the 
behavior of stellar matter. It is also speculated that a transition to a quark phase 
may take place at very high densities, the occurrence of which depends sensitively 
on the properties of the EoS in the hadronic (confined) phase. The presence of 
hyperons in the interior of neutron stars is reported to soften the equation of state, 
with the consequen ce th at the predicted neutron star maximum masses become 
considerably smaller JHQI With recent constraints allowing maximum masses larger 
than previously accepted limits (see previous section), accurate microscopic calcu- 
lations which include strangeness (in addition to other important effects, such as 
those originating from relativity), become especially timely. 

Thus, there are strong motivations for including strange baryons in nuclear mat- 
ter. Moreover, as far as terrestrial nuclear physics is concerned, studies of hyperon 
energies in nuclea r matter natur a lly complement our knowledge of hypernuclei (see, 
for instance, Refs! 111 ' 112 ' 113 1 114 1 ). For example, the EoS of hypermatter is useful 
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in the develop ment of g eneralized mass formulas depending on both density and hy- 
peron fraction !^ ^H^ From the experimental side, additional data are very much 
needed, especially on AA hypernuclei, which would provide information on the A — A 
interaction. Concerning single hypernuclei, analyses of data on A bi nding energies 
constrain the depth of the single-A potential to be 27-30 MevP^ The statu s of 
£ hypernuclei and the E-nucleus potential is more controversial (see 

Ref.™ 

and 

references therein). 



5.2. Hyperons and neutron star matter 

Interacting hypern uclear matter was initially studied within variational 
approach es! ^^ ^ Relat ivistic mean field models were also extended to include 
hyperons! ^0|121 | 122 | 123 | Microscopic mode ls require inclusi on of realistic nucleon- 
hyperon and hypcron-hypcron interactions j 1 24 | 1 25 | 1 26 | 1 27) experimental 
information on these interactions is still scarce. The nucleon-hyperon potentials 
are fitted to AN or EiV scattering data. The information on the h yper on-hyperon 
interaction is limited to the ground state of double-A hypernuclei! I Additional 
constraints can be derived from SU(S) symmetry arguments. 

The common denominator among microscopic calculations of the EoS with hy- 
perons is that they are typically conducted within a non-relativistic framework to- 
gether with r-space local nucleon-hyperon (NY) potentials. Microscopic calculations 
of nuclear matter properties including hyperons have been reported earlier within 
the non-relat ivistic Brueckner-Hartree-Fock framework (BHF) (see, for instance, 
using the Nijmege nucleon-hyperon (NY) meson-exchange po- 
tential. Extensive microscopic work on hyperonic nuclear matter, again within the 
non-relativistic BHF f ramework, has a lso been published by the Barcelona group 
(see, for instance, Refs! 129 ! 130 ! 131 ! 132 !). 

The issue of interacting hyperons and their impact on the 0- s table EoS and 
neut ron star structure has b een confronted, for instance, in 

RefPS] 

using the 

PariPl 

and the Argonne VI jS H nn potentials together with the Nijmegen soft- 
core (NSC89) NY potential! ! 2 ^ A remarkable conclusion from that work is that, 
in the presence of hyperons, the inclusion of nucleonic TBF does not alter the EoS 
appreciably. With nucleons only, TBF bring in considerable repulsion at high den- 
sity which result in a much stiffer EoS. On the other hand, when hyperons are 
present the nucleonic TBF increase the strange baryon population due to the in- 
creased nucleon chemical potential. In turn, this decreases the nucleon population 
with the final net effect on the EoS found to be very small. In other words, it would 
be necessary for the NY interaction to become very repulsive at high density to 
compensate for the loss of nucleons and gain a substantial increase of the star max- 
imum massff^ Results for the neutron star maximum gravitational mass, radius, 
and central density are shown in Table 6. The use of different NN potentials seem 
to produce only small variations. 

More recently, neutron star structure results have been revisited by the Catania 
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Table 6. Neutron star limiting values for different EoS with (Y) and without (no Y) hyperons as 
from Ref. 132 . 



EoS 




/M 


R{k 


m) 


Pc(fm 






no Y 


Y 


no Y 


Y 


no Y 


Y 


V18 


1.64 


1.26 


9.10 


8.70 


1.53 


1.86 


Paris 


1.67 


1.31 


8.90 


8.62 


1.59 


1.84 


V18+TBF 


2.00 


1.22 


10.54 


10.46 


1.11 


1.25 


Paris+TBF 


2.06 


1.26 


10.50 


10.46 


1.10 


1.25 



group with the micr oscopic models we considered earlier in this article ( BHF +TBF 
models from RefP^) along with the NSC89 nucleon- hyperon interaction!!^! There, 
the nucleonic energy densities obtained with the different NN interactions are used 
together with the same NY interaction to calculate stellar matter. It is found that 
the inclusion of hyperons reduces dramatically the maximum mass. It also reduces 
substantially the maximum mass range covered by the models, which goes from 1.8 
to 2.5 Mq with nucleonic energy densities only, see Fig. [18l to the much narrower 
intervals of 1.3-1.4 M Q . Furthermore, the stiffer the original EoS, the larger the 
softening effect from the presence of strangeness. Thus, in this approach, one may 
conclude that hybrid stars (containing hadronic and quark matter) are necessary 
to reproduce the larger mass values consistent with recent observations. 

Overall, there seems to be a consensus among microscopic calculations of the 
BHF type (with or without TBF) that the inclusion of hyperons yields a remark- 
ably low value of a neutron star maximum mass. On the other hand, predictions 
from di fferent ve rsions of relativistic mean field models are quite different from one 
another 

The alternative approach to the EoS with strange baryons presented in the next 
section will suggest that present NY interactions are not yet sufficiently constrained 
to allow for definite conclusions concerning the occurrence of a hybrid phase. 



5.3. A first DBHF calculation of the EoS with nucleons and 
A-hyperons 

It is one purpose of this section to bring in the new aspect of Dirac effects on 
the A hyperon as well as the nucleon. By "Dirac effects" we mean that the single- 
baryon wavefunction is calculated self-consistently with the appropriate effective 
interaction. The origin and nature of these effects on the nucleonic equation of 
state was discussed previously, see Sec. 3.2.1. 

A previous calculation from the Idaho grou of the binding energy of a A 
impurity in nuclear matter showed that Dirac effects on the A hyperon yield a 
moderate reduction of the binding energy. In that calculation, we used the most 
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recent nucleon-hyperon (NY) potential reported in Rrf™ (thereafter referred to 
as NY05) , and observed that G-matrix predictions obtained with NY05 are signifi- 
cantly differerent from cal culat ions using the previous (energy independent) version 

11 OQI 

of the Jiilich NY potentia l I Therefore, in this work we will use both potentials, 
for comparison. As usual, the Bonn B potential's used throughout for the nucleon- 
nucleon (NN) part. 

Previous calculations of the EoS with hyperons have typically been conducted 
within a non-relativistic framework together with r-space local NY potentials. The 
DBHF calculation we describe here uses non-local relativistic momentum-space 
(NN and NY) potentials and a relativistic many-body method, and is therefore 
fundamentally different. Next we give a brief review of the formalism. 

The single-nucleon and single-A potentials are obtained as 

U N (k N ) = U NA {k N ) + U NN (k N ), (63) 

and 

U A (k A ) = U AN (k A ), (64) 

i.e., the AA interaction is neglected. In the equations above, the various terms, 
UbxBz, are the contributions to the potential of baryon B\ from its intcration with 
all baryons of type B2 . They are given by 



U NA (k N )= J2 <f?lufo J+ l\ I ' G NA SJ (k(k N ,k A ),P(k N ,k A ))d 3 k A , 
(Mn + L)(ZSn + l) Jo 



T,L,S,J 



(65) 

k N 

U NN (k N )= ( f^l^?o J+ l\ I " G^ s ' J (k(k N ,k' N ),P(k N ,k' N ))d 3 k' N , 



(66) 



and 



U AN (k A )= £ (2r + 1)(2J + 1} f F G T A ^ J (k(k N ,k A ),P(k Nl k A ))d s k N , 

T.L.S,J ^ A A ' ■ J ° 

(67) 

where the channel isospin T can be or 1 for the NN case and is equal to 1/2 for 
the NA case, and Si,ti (i = N,A) are the spin and isospin of the nucleon or A. 
Notice that 

Una Pa 



(68) 



< U N >= —77^4 / U N {k N )dk 3 N , (69) 



Uan pn ' 

an approximation often used in mean-field approaches. 

The average potential energies of nucleons and A's are determined from 

r" 

Pn{2tt) 3 2 
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1 1 1 



pa 2 j Q 



and 



where the factors of 4 and 2 in Eqs. (|69[) and Eq. (|70[) . respectively, account for 
protons and neutrons in both spin states or A's in both spin states. 
Finally the average potential energy per baryon is obtained as 

< jj > = PN < UN > +PA <UA> , 71 , 
Ptot 

from which, together with a similar expression for the kinetic energy, one obtains 
the average energy per baryon. 

The NA G-matrix is obtained from the Bethe-Goldstone equation 



< NA\G NA (E Q )\NA> = < NA\V\NA> (72) 

+ V <NA\V\NY > — ^— <NY\G N a(E q )\NA>, 

where Eq and E are the starting energy and the energy of the intermediate NY 
state, respectively, and V is an energy-independent NY potential. 

For two particles with masses Mjv and Ma and Fermi momenta kp and kp, 
Pauli blocking requires 

n(ZB\-\ 1 \PP + k\> k F and \aP - k\ > k^ 

^ F >-\ otherwise. (73) 



The above condition implies the restriction 

^ P? + k 2_ [k N ? (^ P) 2 + fc2 _ (fc A )2 

2Pk^ >™sv> , U4j 

where 9 is the angle between the total (P) and the relative (k) momenta of the two 
particles, and M = Ma + M^. Angle-averaging is then applied in the usual way. 

In the present calculation we consider a non-vanishing density of A's but do 
not allow for the presence of real E's in the medium (although both A and E are 
included in the coupled-channel calculation of the NY G-matrix, see Eq. (|73[l). Es- 
sentially we are considering a scenario where a small fraction of nucleons is replaced 
with A's, as could be accomplished by an experiment aimed at producing multi-A 
hypernuclei. Multistrange systems, such as those produced in heavy-ion collisions, 
may of course contain other hyperons on the outset. For small A densities, though, 
as those we consider here, the cascade (S) and the E hyperon are expected to decay 
quickly through the strong processes N + 5 — s- A + A and N + E — > N + A. Under 
these conditions, a mixture of nucleons and A's can be considered "metastable" , 
in the sense of being equilibrated over a time scale which is long relative to strong 
processes. (In fact, the strong reactions mentioned above wo uld h ave to be Pauli 
blocked in order to produce a metastable multistrange system 0301) 
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We neglect the YY' interaction, as very little is known about it. Furthermore, 
non-local momentum-space YY potentials, appropriate for our DBHF framework, 
are not available at this time. For these reasons, we keep the A concentration 
relatively low. 

We have incorporated DBHF effects in the A matter calculation, which amounts 
to involving the A single-particle Dirac wave function in the self-consistent calcu- 
lation through the A effective mass, M£ EE However, a problem that needs to be 
addressed with the Jiilich NY potential in conjunction with DBHF calculations is 
the use of the pseudoscalar coupling for the interactions of pseudoscalar mesons (pi- 
ons and kaons) with nucleons and hyperons. For the reasons described in Sec. 3.2.1 
(that is, the close relationship between Dirac effects and "Z-diagram" contribu- 
tions), this relativistic correction is known to become unreasonably large when 
applied to a vertex involving pseudoscalar coupling. On the other hand, the gradi- 
ent (pseudovector) coupling (also supported by chiral symmetry arguments) largely 
suppresses antiparticle contributions. To resolve this problem, one can make use of 
the on-shell equivalence between the pseudoscalar and the pseudovector coupling, 
which amounts to relating the coupling constants as follows: 

9ps fpv 7 (^) 

Trip S 

where g ps denotes the pseudoscalar coupling constant and f pv the pseudovector 
one; m ps , Mi, and Mj are the masses of the pseudoscalar meson and the two 
baryons involved in the vertex, respectively. This procedure can be made plausible 
by writing down the appropriate one-boson-exchange amplitudes and observing 
that, redefining the coupling constants as above, we have (see RefPJ for the two- 
nucleon case) 

V pv = V ps + (76) 

where the ellipsis stands for off-shell contributions. Thus, the pseudoscalar coupling 
can be interpreted as pseudovector coupling where the off-shell terms are ignored. 
This is what we apply in our DBHF calculations. More concretely, from a given 
pseudoscalar potential, first we extract the corresponding f pv from Eq. ([75)1 . In 

the medium, V* s will contain the coupling constant g = — '— — -f pv - When the 
baryon masses are reduced through the effective mass prescription, the reduction 
of the masses appearing at the denominator of the momentum-dependent part of 
the potential will be "balanced" be the equally reduced masses at the numerator of 
the coupling constant. This prescription prevents the OBE diagram from growing 
to unreasonably large values when effective masses are employed. 

Next, we will be showing results for both the NY05 potential I 138 | and the pre- 
vious versi on of the Jiilich NY potential, NY94, specifically the model referred to 
as A in RefPSl. 

In Fig. [23] we show the energy per particle as a function of density for different A 
concentrations as obtained from DBHF calculations along with the NY94 potential. 
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As more nucleons are replaced with A's, generally less binding energy per particle 
is generated. This is due to the weaker nature of the NA interaction relative to 
the NN one. Furthermore, the A Fermi momentum grows rather quickly with p\, 
since only two A's can occupy each state, rather than four, which implies a fast rise 
of the hyperon kinetic energy. (Although, for small hyperon densities, there is at 
first some reduction of the kinetic energy, due to the fact that A's have larger mass 
and can occupy lower energy levels.) We notice that the saturation density remains 
essentially unchanged with increasing hyperon concentrations. As density grows, 
though, larger A concentrations start to yield increased attraction. One must keep 
in mind that, especially at the higher densities, the NN interaction become less and 
less attractive due to medium effects, in particular repulsive Dirac effects. Thus, 
removing nucleons from the system can actually amount to increased binding. 

Moving now to Fig. [24l where the NY05 potential is adopted, we see a very 
different scenario. We recall that the NY05 mode l is cons iderably more attractive, 

11 ^711 49I 

yielding about 50 MeV for the A b inding energy 1 J ' 1 1 whereas a value close to 
30 MeV was found with NY94 11321. Naturally, we expect these differences to reflect 
onto the respective EoS predictions. Here, increasing the hyperon population yields 
more binding, a trend opposite to the one seen in the previous figure. Again, the final 
balance is the result of a combination of effects. The fact that the NN component 
is repulsive at the higher densities, together with the more attractive nature of 
the NY05 potential, determines here a net increase in attraction with decreasing 
nucleon density. However, the additional binding becomes smaller and smaller with 
increasing hyperon concentration, indicating that, at sufficiently large A densities, 
the net balance may turn repulsive. Notice also that in Fig. [M] the minimum moves 
towards higher densities, signifying that baryon pairs favor a smaller interparticle 
distance. 

To summarize, one must keep in mind that the energy/particle is the result of 
a delicate balance of both the kinetic energies and the potential energies of the two 
baryon species, being weighed by the respective densities. Thus, although the NN 
interaction is generally more attractive than the NA one, the net effect of replacing 
nucleons with A's will depend sensitively on the nature of the NY interaction that's 
being put into the system, as well as the "stiffness" of the original, nucleonic, EoS. 

This is further confirmed in Fig. 1251 where we show predictions from conven- 
tional BHF calculations (i.e., no "Dirac" effects). As in the previous figures, the 
red curve is the nucleonic EoS. The solid(dashed) curves are the predictions with 
NY94(NY05). Clearly, the effects are opposite depending on the NY interaction. 
Qualitatively, the trend seen in each group of curves (solid or dash) is approximately 
consistent with the one observed previously in the corresponding (i. e. same NY po- 
tential) DBHF predictions. However, the "cross over" of the curves seen in Fig. [23] 
at about twice saturation density is due to a large extent to the more repulsive 
nature of the nucleonic EoS in the DBHF calculation (see comments above). 

In conclusion, the effect of hyperons on the EoS is strongly dependent upon the 
baseline (nucleonic) EoS as well as the NY potential model. With regard to the 
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Fig. 23. (Color online) Energy/particle as a function of density in symmetric nuclear matter 
for various A concentrations Yj±. Predictions obtained from DBHF calculations with the NY94 
potential. 



first issue, we stress the importance of starting from a realistic EoS, such as the one 
predicted in DBHF calculations or obtained with the incl usion of three-body forces. 
The second observation confirms the conclusions of 

RefPSl 

concerning the large 

uncertainties originating from the bare NY potentials. Unfortunately, the existing 
data do not set sufficient constrains on the potentials, as demonstrated by the fact 
that different potentials can fit the availa ble s cattering data equally accurately but 
produce very different scattering lengths 0421 

Of course, the actual fraction of hyperons present in star matter must be deter- 
mined by the equations of /3 stability and charge neutrality for highly asymmetric 
matter containing neutrons, protons, hyperons, and leptons as described in Sec. 5.1. 
What we have learnt at this time is that the predicted energy/particle in symmetric 
matter is very sensitive to the chosen NY interaction. The uncertainties due to the 
model dependence discussed in this paper are likely to impact any conclusions on 
the properties of strange neutron stars, which therefore must be interpreted with 
caution. These include considerations of deconfmement and possible transition from 
hadronic to quark matter, which depend sensitively on the equation of state in the 
hadronic phase. 



5.4. Neutron star matter with non-interacting hyperons 

As anticipated at the end of Sec. 5.2, the uncertainties that we have demonstrated in 
Sec. 5.3 suggest that present NY interactions are not sufficiently constrained to al- 
low for robust conclusions concerning the occurrence of a hybrid phase. Therefore, to 
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Fig. 24. (Color online) Energy/particle as a function of density in symmetric nuclear matter 
for various A concentrations Predictions obtained from DBHF calculations with the NY05 
potential. 





Baryon density (fm ) 



Fig. 25. (Color online) Energy /particle as a function of density in symmetric nuclear matter for 
various A concentrations. Predictions obtained from BHF calculations. The solid(dashed) lines are 
obtained with the NY94(NY05) potentials. 



avoid those uncertainties, we will now try to estimate the impact of hyperons on the 
/3-stable DBHF EoS considering non-interacting hyperons. Equations (|59I60I62|) are 
solved (for fixed total baryon density) treating A and E~ as free (non-relativistic) 
fermions. The nucleonic energies are taken from our DBHF EoS for IANM, eval- 
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p (fnr 3 ) 

Fig. 26. (color online) Various particle fractions in /3-stable matter. Neutron: long-dashed black; 
proton: solid black; S — : dash-dotted green; A: dash-double-dotted blue; electrons: short-dashed 
red; muons: dotted purple. 

uated at some nucleoli density, p n + p v , with p n and p p to be determined along 
with the other particle fractions. (To allow for the solution of the 6x6 (including 
the normalization condition) system of algebraic equations for the six unknown 
particle concentrations, our nucleonic EoS is first parametrized in terms of analytic 
functions of density.) T hose fractions are shown in Fig. 1261 The trend is very similar 
to the one shown in with free hyperons, which, in turn , is q ualitatively con- 

sistent with the one seen when the NY interaction is turned on|122l The onset of the 
E~ baryon occurs at a density of about 0.3/m -3 . Leptons start to disappear from 
the system after hyperon formation, an indication that charge neutrality through 
£~ production is energetically more favorable than through j3 decay. The onset of 
A production occurs at a density of approximately 0.6/m -3 . These thresholds are 

II QoJ 

in good agreement with those reported in Ref. for both free and interacting hy- 
perons, when the BHF+TBF predictions are considered (the TBF shifts down the 
threshold densities for A and E~ formation as compared with the case of two-body 
forces only). 

In Fig. [27l we show the pressure in /3-stable matter with the Idaho DBHF EoS 
with and without (fre e) hy perons. The differences are dramatic but consistent with 

[Tool 

previous observations.^^ The effect of including the NY interaction is quite small 
compared with the differences seen in Fig. 1271 Thus, one conclusion seems appro- 
priate at this point: Regardless the chosen NY interaction, the effect of adding 
strange baryons is a dramatic softening of the EoS, (consistent with a large re- 
duction of the neutron star maximum mass) . Furthermore, given the uncertainties 
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Fig. 27. (color online) Pressure in /3-stable matter from the DBHF calculation with no hyperons 
(black dash) and with free hyperons (solid red). 



arising from the NY potential demonstrated in the previous section, the effect of 
free hyperons may very well provide a realistic, average estimate of the impact of 
including strangeness. 

In Fig. [28] we compare energy densities in /3-stable matter as obtained from 
the Idaho DBHF EoS wit h fre e hyperons and from the two models (AVI 8 + TBF, 
Paris + TBF) used in RefffSSl along with the NSC NY potential!^ The similarity 
between the curves indicates that the effect of including the NY interaction in star 
matter is relatively small and/or that differences in the respective nucleonic EoS 
may compensate those which originate from the treatment of hyperons. 

The neutron star maximum mass we obtain with the DBHF EoS and free hy- 
perons is close to 1.2 M Q . Considering that inclusion of the NY interaction would 
increase the pressure at high density, this result is not inconsistent with those shown 
in Table 6. In summary, our analysis confirms a large softening of the /3-stable EoS 
as a consequence of including strange baryons, accompanied by a strong reduction 
of the maximum mass, a conclusion that appears to be nearly model independent. 
The softening is mostly caused by conversion of the kinetic energy of the existing 
particle species into masses of the newly produced particles. The resulting maxi- 
mum mass reduction appears in contraddiction with present constraints, and thus 
leaves intriguing questions open concerning phase transitions and degrees of free- 
dom appropriate for high-density stellar matter. 
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Fig. 28. (color online) Energy density in /3-stable matter from various models: The black dashed 
and dotted lines (essentially identical) are the results of BHF calculations with the Argonne V18 
and the Paris potentials, respectively, including TBF and hyperons interacting via the Nijmegen 
soft-core NY potential (see Ref. 132 ). The red solid line is obtained with the DBHF calculation for 
the nucleonic EoS and free hyperons. 

6. Conclusions and Outlook 

The EoS of hadronic matter enters in a variety of systems, from atomic nuclei to 
the most compact form of matter found in the universe, namely the matter in the 
interior of neutron stars. EoS-related properties can and should be calculated ab 
initio and consistently from the same basic nuclear forces. They include: the nuclear 
matter "optical potential", effective in-medium cross sections, nucleon mean free 
path in nuclear matter, neutron skins, and neutron star properties. 

We have reviewed recent calculations of the EoS, over a large range of densities 
and isospin asymmetries, with a particular eye on the microscopic approach. The 
obvious advantage of the latter lies in the opportunity of interpreting the predic- 
tions in terms of the input nuclear forces and their behavior in the many-body 
environment. In turn, this facilitates the identification of potentially missing dy- 
namics. 

We have discussed differences among microscopic models which do or do not 
include explicit TBF. It is quite clear that model dependence amongst predictions 
can be quite large. Naturally, this is especially the case at those densities where 
constraints are the weakest. 

Rich and diverse effort is presently going on to improve the available constraints 
on the EoS or find new ones. These constraints are usually extracted through the 
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analysis of selected heavy-ion collision observables. At the same time, partnership 
between nuclear physics and astrophysics is becoming increasingly important to- 
wards advancing our understanding of exotic matter. Perhaps the best prospects to 
discriminate among families of EoS, especially at high density, are in more accurate 
measurements of neutron stars radii and/or moments of inertia. 

Other EoS-related issues of contemporary interest which we have not reviewed 
here include temperature dependence and polarizability of nuclear/neutron matter. 
Knowledge of the finite-temperature EoS plays a crucial role in the final stages 
of a supernova evolution. Microscopic models which can successfully describe the 
ground state of nuclear matter should be the starting point to move on to its excited 
states. 

As a final note, the emergence of FRIB 21 will mark a turning point in the 
progress of experimental nuclear physics, possibly unveiling unknown areas on the 
nuclear chart. It is thus imperative that such large-scale projects be constantly 
supported by theoretical calculations with predictive power. 
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